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ABSTRACT 
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I.  INTRODUCTION 

A.  BACKGROUND 

This  thesis  was  a  continuation  of  the  work  currently  being  performed  at  the 
Naval  Postgraduate  School  and  was  supported  by  the  contributions  attained  in  Ref- 
erences 1,  2,  3  and  4.  The  goal  of  this  research  was  to  obtain  data  concerning  the 
effect  of  propellant  properties  upon  the  performance  of  soHd  propellant  rocket  mo- 
tors. Data  analysis  is  performed  on  holographic  images  obtained  from  combustion 
products  created  during  rocket  firings.  Holograms  are  obtained  through  the  use  of 
a  pulsed  ruby  laser  and  a  holocamera.  Reconstruction  of  the  holographic  image  for 
analysis  is  achieved  with  a  krypton  laser,  variable-power  microscope,  and  video  cam- 
era. Unfortunately,  the  optical  diffuser  required  to  record  the  hologram  creates  a 
speckle-corrupted  image  which  must  be  filtered  before  the  data  retrieval  process  is 
initiated.  The  data  attained  from  these  procedures  is  used  in  performance  studies, 
stability  analysis  predictions,  and  exhaust  signature  diagnosis. 

B.  THESIS  OBJECTIVES 

The  thesis  objectives  were: 

1.  To  implement  software  modifications  in  existing  local  code  in  support  of  the 
new  computer  system  configuration. 

2.  To  enhance  program  code  portability. 

3.  To  optimize  program  code  and  attempt  to  reduce  user  interaction  where  possi- 
ble. 

4.  To  create  a  menu  format  for  executable  files.  ■  ' 

All  objectives  were  accomplished.  ITEX/PCp/u5  library  routines  have  been  suc- 
cessfully integrated  into  the  image  analysis  programs.  The  program  file  genfunc.c  was 
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used  to  optimize  program  size,  and  new  functions  were  created  for  user  convenience. 
A  menu  format  was  devised  to  make  the  image  analysis  procedure  easier  to  use. 

C.     SOFTWARE  SUPPORT 

The  associated  system  software  consists  of  the  following  items: 

•  Image  analysis  programs  created  via  thesis  research  in  References  1,  2,  3,  and 

4, 


Imaging  Technology's  ImageActionp/u5  Software  package, 
Imaging  Technology's  ITEX/PCp/u5  software  library, 
IBM  Disk  Operating  System  version  3.3  (DOS), 

•  Microsoft  C  optimizing  compiler  version  5.0  (MSC),  and 

•  Microsoft  Codeview  version  2.10. 

The  ImageAction;?/us  software  package  is  a  menu-driven  system  that  allows  the 
user  to  process  stored  or  Hve  video  images  [Ref.  5].  The  two-level  menu  organization 
provides  image  analysis,  graphics,  geometric  operations,  and  image  processing  capa- 
bihties.  This  software  package  also  has  the  ability  to  utilize  user-generated  special 
purpose  script  files.  The  script  files  are  created  in  C  language  and  are  supported  by 
libraries  created  by  Imaging  Technology,  Inc.  Imaging  Technology's  ITEX/PCp/u5 
software  hbraries  are  compatible  only  with  the  Microsoft  PASCAL,  FORTRAN,  and 
C  compilers.  These  hbraries  consist  of  image  processing  and  graphics  function  sub- 
routines used  in  conjunction  with  the  PCVISI0Np/u5  Frame  Grabber  board,  a  frame 
memory  and  video  digitizer  installed  in  the  IBM  PC/AT. 

D.     THESIS  ORGANIZATION 

This  thesis  is  organized  into  five  chapters.  The  introductory  chapter  provides 
the  project  background  and  goals.  It  also  discusses  thesis  objectives  and  organization. 


Chapter  II  describes  the  system  hardware  modifications.  The  PCVISION  and 
the  PCVISIONp/t/5  hardware  configurations  are  compared,  and  the  respective  hard- 
ware capabilities  are  presented.  Also  the  hardware  performance  is  briefly  discussed. 

The  third  chapter  gives  a  synopsis  of  the  image  processing  procedure  and  dis- 
cusses actual  software  programs  and  improvements  at  length.  The  chapter  also  de- 
scribes the  techniques  and  software  development  tools  utihzed  for  program  modifica- 
tion. Finally,  Chapter  III  provides  program  specifics. 

The  software  performance  analysis  is  provided  in  Chapter  IV.  Program  analysis 
consisted  of  comparisons  in  the  following  areas:  executable  program  size,  program 
execution  times,  and  program  result  verification. 

Lastly,  Chapter  V  summarizes  the  findings  obtained  during  the  course  of  this 
research  and  provides  recommendations  for  future  thesis  research. 


II.  SYSTEM  HARDWARE  MODIFICATIONS 

A.     HARDWARE  CONFIGURATION 

Digital  image  processing  involves  the  achievement  of  the  following  objectives 
through  a  computer  interface:  recognition,  segmentation,  enhancement,  or  analysis 
of  an  image.  The  image  processing  and  digitization  system  consists  of  the  following 
components: 

•  An  IBM  PC/AT  computer  with  40  megabyte  hard  drive, 

•  An  Intel  INBOARD  386/AT  microprocessor  and  80387  math  coprocessor, 

•  An  AST  ADVANTAGE  memory  expansion  board, 

•  A  PCVISI0Np/w5  Frame  Grabber  board, 

•  A  computer  monitor  with  EGA  display, 

•  A  Panasonic  color  video  data  monitor,  and 

•  A  video  cassette  recorder. 

A  block  diagram  of  the  basic  elements  of  the  image  processing  system  is  presented 
in  Figure  2.1.  One  of  the  primary  components  of  this  image  processing  system  is 
the  video  cassette  recorder,  which  provides  the  introduction  of  the  image  into  the 
PCVISI0Np/w5  Frame  Grabber  board.  This  image  signal  consists  of  timing  and 
analog  video  information  in  a  standard  RSI 70  video  format.  Through  the  digitization 
process,  the  analog  signal  is  converted  to  a  digital  format  for  processing  and  storage. 
The  IBM  PC/AT  central  processing  unit  accesses  these  digital  values  via  the  frame 
memory. 

The  PCVISI0Np/u5  Frame  Grabber  is  a  single  board  that  is  connected  directly 
into  one  of  the  16-bit  expansion  slots  of  the  IBM  PC/AT  computer.  Pixels  are  stored 
in  frame  memory  after  the  Frame  Grabber  digitizes  the  analog  video  signal  at  a 


1 


>- 

UJCC 

2o 

<5 

^2 

/ 

\ 

i 

o 

o 

a: 

o 

LU 

_l 

N 

Q 

> 
< 

_l 

Q. 
CO 

Q 

1 

/ 

\ 

\ 

/ 

< 

Jr    . 

_J 

<cr 

lU 

zO 

S 

cc^- 

< 

UJ^ 

O 

^1 

y   V 

/     \   1 

J 

Figure  2.1:  Elements  of  the  Image  Processing  System.   The  Dotted  Rect- 
angle Encloses  Frame  Grabber  Board.   (From  Ref.  7) 
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Figure  2.2:  Display  of  AVW  and  Hidden  Areas  of  Dual  Memory  Configu- 
ration.  (From  Ref.  6) 

rate  of  thirty  frames  per  second  with  eight  bits  of  accuracy.  In  previous  work,  a 
PCVISION  Frame  Grabber  board  was  utilized  (versus  the  PCVISlONp/us  version). 
The  PCVISI0Np/u5  board  permits  either  two  512  x  512  pixel  images  or  one  640  x 
512  pixel  image  to  be  stored  within  frame  memory  [Ref  7:p.  1-4].  In  the  PCVISION 
board,  only  one  512  x  512  pixel  image  can  be  stored  [Ref.  8:p.  1-2].  The  bottom 
thirty-two  lines  of  the  video  output  from  both  versions  of  the  Frame  Grabber  board 
are  within  a  "hidden  area"  which  can  be  accessed  by  the  host  computer  but  not  viewed 
on  the  video  monitor  without  scrolling  the  image.  The  active  video  window  (AVW) 
is  the  area  within  frame  memory  in  which  the  pixel  values  are  displayed  or  stored  via 
the  scanning  circuit  [Ref  6:  p.  1-14].  Figure  2.2  and  Figure  2.3  indicate  respectively 
the  dual  and  single  memory  configuration  hidden  and  AVW  areas.  Currently,  only 
the  dual  memory  configuration  is  utilized  in  thesis  research. 
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Figure  2.3:  Display  of  AVW  and  Hidden  Areas  of  Single  Memory  Config- 
uration.  (From  Ref.  6) 

Prior  to  display  of  the  image  on  the  video  monitor,  look-up  tables  (LUT)  are 
employed  for  the  transformation  of  pixel  intensities.  The  original  PCVISION  Frame 
Grabber  had  only  four  LUTs  associated  with  output  [Ref.  8:p.  2-8].  The  current 
Frame  Grabber  board  has  one  input  digitization  path  and  three  output  channels,  each 
allocated  eight  LUTs  [Ref.  7:p.  1-4].  These  LUTs  allow  for  real  time  processing  by 
virtue  of  their  ability  to  afford  simple  point  transformations  of  the  256  gray  levels 
without  any  processing  delay.  Look  up  tables  are  addressed  by  eight  bits  of  pixel 
information  provided  by  the  frame  memory.  A  digital-to-analog  converter  and  the 
look  up  tables  constitute  the  display  logic  unit  of  both  versions  of  the  Frame  Grabber 
boards.  The  digital-to-analog  converter  receives  gray  intensity  values  corresponding 
to  the  actual  value  addressed  within  the  LUTs  [Ref.  8:p.  2-9]. 


Both  versions  of  the  Frame  Grabber  board  are  compatible  with  the  IBM  PC/AT 
computer  bus  structure.  This  host  interface  provides  the  communication  path  be- 
tween the  Frame  Grabber  boards  and  the  personal  computer.  The  previous  PC- 
VISION  board  host  interface  consisted  of  three  components:  the  Interrupt  Logic, 
Control  Register  Interface,  and  the  Frame  Memory  Interface  [Ref.  8:p.  2-1].  The 
PCVISI0Np/u5  host  interface  is  composed  of  the  Control  Register  Interface  and  the 
Frame  Memory  Interface  [Ref.  7:p.  4-2].  The  Interrupt  Logic  had  afforded  the 
PCVISION  Frame  Grabber  the  ability  to  be  controlled  in  real  time  while  the  IBM 
personal  computer  processed  other  instructions  [Ref.  8:p.  2-3].  This  outdated  design 
feature  has  been  superceded  by  transparent  computer  access  to  frame  memory  in  the 
PCVISI0N>/?Z5  system  [Ref.  7:p.  1-6]. 

Within  the  input/output  (I/O)  space  of  the  PC/AT,  sixteen  bytes  are  reserved 
for  the  PCVISI0Np/u5  Frame  Grabber  control  registers.  Currently,  only  twelve  of 
these  sixteen  bytes  are  being  used  by  the  system.  Access  to  the  remaining  bytes 
is  prohibited;  they  are  reserved  for  future  system  growth  [Ref.  7:p.  4-3].  On  the 
contrary,  the  PCVISION  board  requires  32  bytes  within  the  I/O  space  for  its  control 
registers.  Only  the  first  seven  of  these  bytes  are  actually  utilized  for  control  functions 
[Ref.  8:p.  2-3]. 

The  frame  memory  of  the  PC VISIONp/«s  system  is  apportioned  into  eight  equiv- 
alent blocks  of  64K  bytes.  Block  dimensions  are  512  pixels  by  128  pixels,  and  the 
physical  layout  of  the  dual-store  configuration  frame  memory  is  shown  in  Figure  2.4. 
This  configuration  can  be  contrasted  to  that  of  the  640x512  image  in  which  memory 
is  apportioned  into  eight  blocks  of  1024  pixels  by  64  pixels.  The  single-store  frame 
memory  configuration  can  be  viewed  in  Figure  2.5.  Individual  selection  of  these  blocks 
is  controlled  by  select  bits  within  the  Control  Register.  As  viewed  by  the  host,  each 
block  gives  the  semblance  consecutive  data  bytes  originating  at  the  memory  base 
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Figure  2.4:  PCVISIONp/u5  Dual  Frame  Memory  Layout.   (From  Ref.   7) 

address  [Ref.  7:p.4-3].  For  the  PC  VISION  system,  the  frame  memory  was  subdivided 
into  four  blocks  of  256  pixels  by  256  pixels,  forming  a  quadrant  structure  as  displayed 
in  Figure  2.6.  At  any  time  only  one  block  can  be  accessed  from  the  personal  computer 
bus  [Ref.  8:p.  2-7]. 

Both  versions  of  the  Frame  Grabber  board  require  the  assignment  of  a  single 
block  of  64K  bytes  within  the  IBM  PC/AT  memory  space  located  above  640K.  Spe- 
cialized software  controls  the  mapping  of  individual  blocks  of  frame  memory  into  this 
reserved  space  [Refs.  7  and  8].  In  the  case  of  the  PCVISION;7/u5  board,  the  memory 
base  address  is  factory-configured  at  address  AOOOOII.  This  address  may  be  modified 
upon  board  installation  provided  the  subject  address  selected  is  within  the  guidelines 
set  forth  in  Reference  7. 

It  was  originally  intended  that  this  thesis  work  would  be  implemented  on  a 
COMPAQ  DESKPRO  386/20  computer  system.   The  PCVISI0N;?/n5  memory  base 
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address  was  configured  to  an  approved  value  of  DOOOOH.  Therefore,  the  address  region 
from  DOOOOH  to  DFFFFH  was  reserved  within  the  computer  memory  space,  and  all 
other  memory  and  peripherals  were  restricted  access  to  that  memory  region.  Upon 
subsequent  testing  and  diagnostic  checks,  it  was  found  that  there  was  an  unresolved 
memory  conflict  within  the  COMPAQ  computer  memory  which  affected  some  of  the 
operations  of  the  Frame  Grabber  board.  Due  to  necessity,  the  board  was  consequently 
installed  in  the  IBM  PC/AT.  Memory  base  address  DOOOOH  was  again  selected;  this 
time  there  were  no  resident  memory  conflicts. 

B.     HARDWARE  PERFORMANCE 

With  the  introduction  of  the  PCVISIONp/us  Frame  Grabber  board,  the  abihty 
to  process  two  separate  images  in  frame  memory  is  invoked.  The  dual  memory  con- 
figuration allows  the  host  to  select  either  image  area  for  analysis  and  manipulation. 
A  display  memory  function  allows  the  user  to  view  independently  either  image  area. 
These  capabilities  allow  for  considerable  future  flexibility  in  the  analysis  of  the  rocket 
motor  holograms. 

The  PCVISI0Np/w5  board  not  only  provides  more  image  processing  capabilities 
than  the  PCVISION  board,  but  it  also  establishes  the  need  for  additional  function 
libraries  to  support  its  abilities.  For  this  reason,  image  analysis  programs  devised  in 
prior  thesis  research  had  to  be  modified  to  reflect  the  new  supporting  Hbraries.  These 
programs  and  their  resulting  modifications  are  emphasized  in  the  following  chapter. 


11 


III.  SOFTWARE  PROGRAMS  AND 
IMPROVEMENTS 

A.  IMAGE  PROCESSING  PROCEDURE 

As  was  discussed  in  the  introductory  chapter,  holographic  techniques  are  em- 
ployed to  obtain  the  image  of  combustion  products  from  small  rocket  motors  during 
firing.  During  the  hologram  reconstruction  process,  a  krypton  laser  is  used  in  con- 
junction with  a  variable  power  microscope.  Due  to  the  presence  of  a  difFuser  that  is 
required  to  ehminate  thermal  gradient  effects,  speckle  noise  is  introduced  into  the  cap- 
tured image.  The  following  paragraphs  describe  the  procedure  developed  to  acquire 
and  process  this  image  for  statistical  study. 

1.     Image  Acquisition 

The  speckle-corrupted  image  is  recorded  for  future  processing  via  a  combi- 
nation of  a  video  cassette  recorder  and  low-light-level  camera  attached  to  the  variable 
power  microscope.  The  magnification  power  of  the  microscope  determines  the  rel- 
ative size  of  the  pixel  elements  comprising  the  image.  Pixels  are  the  standard  unit 
for  the  measurement  of  the  features  within  the  image.  Object  sizes  are  determined 
by  scaling  the  measured  object  width  (in  pixels)  by  the  dimension  spanned  by  one 
pixel  (in  microns  per  pixel).  A  quantization  error  is  introduced  into  this  calculation 
due  to  the  necessity  of  measuring  dimensions  with  an  integer  number  of  pixels.  As 
the  width  decreases,  the  proportional  level  of  uncertainty  increases,  which  produces 
a  more  erroneous  scale  factor  calculation.  A  calibration  object  is  inserted  within  the 
holographic  image  to  determine  the  width  spanned  by  a  pixel.  It  is  essential  that  this 
object  covers  an  appreciable  portion  of  the  monitor  screen  (at  least  one  half  of  the 
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area)  to  counteract  this  quantization  error.  The  calculation  of  the  conversion  factor 
from  pixels  to  microns  becomes  increasingly  important  in  subsequent  steps  of  image 
processing. 

2.  Image  Digitization 

From  the  VCR  tape  obtained  from  the  hologram  reconstruction,  the  menu- 
driven  ImageActionp/us  software  is  used  to  "grab"  and  digitize  the  desired  image. 
The  PCVISI0N/>/u5  Frame  Grabber  board  implements  this  procedure.  The  digitized 
video  image  is  comprised  of  a  512  x  480  array  of  pixels  having  gray  intensity  levels 
ranging  from  0  to  255.  These  intensity  levels  extend  from  blackest- black,  level  0,  to 
whitest- white,  level  255. 

3.  Speckle  Reduction  Filtering 

Image  filtering  is  the  paramount  technique  used  in  speckle  reduction.  The 
amount  of  speckle  reduction  achieved  by  filtering  algorithms  is  gauged  by  a  calculation 
of  the  speckle  index.  Speckle  index  is  defined  as  the  ratio  of  the  image's  local  deviation 
to  the  local  mean  of  pixel  values,  and  this  measure  is  indicative  of  the  random  speckle 
noise  inherent  in  the  image  [Ref.  l:pp.  17-19]. 

The  three  filtering  algorithms  are  based  upon  synthetic  aperture  radar 
(SAR)  theory  presented  in  References  9  and  10.  There  are  three  separate  filter  mod- 
els. Derived  from  nonlinear  techniques,  two  filter  designs  are  based  upon  statistical 
methods.  The  third  filter  model  has  its  foundation  in  a  geometric  hulling  algorithm. 
The  geometric  hulling  algorithm  produces  the  best  results  by  enabling  the  user  to 
discontinue  filtering  at  a  less  severe  level  of  resolution  degradation. 

4.  Image  Threshold 

In  an  ideal  image,  feature  particles  and  the  background  would  have  distinct 
gray  intensity  levels.    The  process  of  thresholding  creates  a  binary  image  in  which 
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image  features  will  appear  black  and  the  background  white.  In  reality,  the  thresholded 
image  will  contain  some  amount  of  speckle  manifested  &s  feature  particles  and  will 
have  lost  portions  of  actual  particles  during  the  filtering  and  thresholding  process. 
Image  evaluation  errors  are  inherent  whenever  the  grey  levels  of  the  object  overlap 
the  gray  levels  of  the  image  background.  The  selection  of  a  threshold  value  is  obtained 
from  observation  of  the  filtered  image.  ITEX/PCp/u5  software  allows  the  user  to  select 
a  threshold  value  for  which  the  designated  LUT  is  modified.  This  modification  is  not 
made  permanent  until  the  image  is  mapped  into  frame  memory,  allowing  the  user 
to  rapidly  iterate  values  to  make  a  determination  on  the  most  appropriate  threshold 
value. 

5.  Feature  Identification 

Once  the  filtered  image  has  been  thresholded,  it  is  ready  to  undergo  fea- 
ture identification.  It  is  imperative  that  the  programmed  image  has  been  previously 
thresholded  due  to  the  logic  of  the  identification  algorithm.  This  procedure  locates 
all  connected  adjacent  feature  pixels  and  assigns  feature  value  numbers  to  them. 
Through  this  labeling  process,  interconnected  pixels  are  identified  as  an  individual 
object.  The  total  number  of  objects  recognized  serves  as  input  for  the  next  step 
of  image  processing.  This  number  is  restricted  to  32,767  or  below  based  upon  the 
limitations  imposed  by  the  C  language  on  integer  values  [Ref.  4:p.  9]. 

6.  Feature  Sizing 

The  procedure  of  feature  sizing  determines  both  the  maximum  vertical  and 
horizontal  diameters  of  the  identified  objects.  The  total  area  of  the  subject  particles 
is  computed,  and  all  these  measurements  are  delineated  in  terms  of  microns.  As 
described  before,  a  conversion  factor  is  necessary  for  manipulation  of  the  feature 
particles  pixel  measurements  into  actual  physical  dimensions.  This  constant,  called 
SCALE-FACTOR,  is  presently  based  upon  the  dimensions  of  a  threaded  screw  at  the 
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same  magnification  of  the  captured  image.  The  screw  threads  serve  as  a  calibration 
object  to  transform  image  object  pixel  size  data  into  actual  physical  measurement 
criteria.  The  final  output  of  the  feature  sizing  routine  is  presented  in  tabular  form 
suitable  for  input  into  the  STATGRAPHICS  software  package,  a  statistical  analysis 
tool. 

7.      Histogram  production 

To  exhibit  feature  size  data  in  a  histogram  form  for  ready  evaluation,  the 
STATGRAPHICS  software  package  is  utilized.  This  particular  software  package  was 
selected  because  it  is  compatible  with  the  IBM-PC/AT  system  and  is  relatively  easy 
to  use.  Histograms  produced  by  STATGRAPHICS  are  used  for  further  statistical 
analysis  to  determine  total  particle  distributions  of  holograms. 

B.     PROGRAM  OPTIMIZATION 

The  primary  purpose  of  this  thesis  research  was  to  modify  existing  C  language 
routines  to  make  them  compatible  with  the  ITEX/P Cplus  software.  Subsequent  pro- 
gram optimization  was  performed  in  support  of  this  goal.  One  significant  improve- 
ment was  the  removal  of  all  general-use  functions  from  the  main  programs.  Although 
a  general  support  functions  program  file  had  been  created  in  the  previous  rendition, 
this  file  had  not  been  adequately  utilized.  As  a  result,  all  the  main  programs  had 
needlessly  large  file  sizes. 

The  C  language  provides  the  ability  to  use  the  infinitely-abusable  GO  TO  state- 
ment. This  creates  a  situation  in  which  branching  occurs  in  a  deeply  nested  structure. 
The  GO  TO  organization  may  prove  effective  in  some  instances,  but  it  always  pro- 
duces confusion  for  the  individual  having  the  misfortune  to  try  to  modify  that  code 
sequence.  For  this  reason,  program  structures  were  modified  to  eliminate  all  GO  TO 
branch  statements,  thereby  making  code  more  modular  and  portable. 
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A  menu  format  was  developed  for  use  with  the  executable  programs.  This 
format  will  facihtate  demonstrations  of  image  processing  and  makes  the  routines 
"user- friendly" .  In  addition,  supporting  functions  were  devised  for  some  of  the  pro- 
grams to  further  ehminate  the  need  for  user  interaction  with  the  programs  once  the 
selected  image  had  been  loaded  and  image  processing  regions  were  identified.  These 
modifications  facihtated  more  rapid  program  execution.  A  more  detailed  discussion 
of  program  modification  specifics  will  occur  in  a  succeeding  section. 

C.     SOFTWARE  DEVELOPMENT  TOOLS 

The  Microsoft  Codeview  Debugger  program  was  utihzed  for  program  develop- 
ment. It  requires  that  specific  compiler  directives  be  used  during  program  compi- 
lation. Although  Codeview  is  an  excellent  developmental  tool,  program  debugging 
can  become  a  tedious  process.  No  program  modifications  can  be  performed  while  in 
Codeview,  and  for  that  reason,  it  is  not  an  adequate  tool  for  what  is  needed  for  pro- 
gram modules  of  this  size.  The  ITEX/PCp/ws  routines  require  large  memory  module 
handlers,  and  few  other  debuggers  are  able  to  meet  this  requirement.  Quick  C  version 
2.0  is  capable  of  handling  large  memory  modules  and  has  an  online  source  debugger 
which  is  more  versatile  than  Codeview  [Ref.  11].  It  is  recommended  that  in  future 
work  with  these  program  modules,  Quick  C  should  be  included  in  system  upgrades. 

Several  batch  files  were  created  for  compilation  and  linker  processes  since  Code- 
view  required  special  compiler  directives.  There  are  two  different  versions  of  these 
files:  compilation  and  link  with  debugger  information,  and  compilation  without  de- 
bugger information  for  file  optimization  and  compaction.  Listings  and  justifications 
for  these  DOS  Batch  files  are  referenced  in  Appendix  K. 
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D.     PROGRAM  SPECIFICS 

The  succeeding  paragraphs  discuss  the  programs  that  were  developed  in  support 
of  this  thesis  research.  These  existing  programs  were  modified  to  be  compatible  with 
the  ITEX/FCplus  software  hbraries  and  the  ImageActionp/us  routines.  Where  at 
all  possible,  program  optimization  was  performed.  In  some  ca^es,  code  had  to  be 
completely  revamped  to  accommodate  new  requirements.  All  source  code  hstings  are 
included  as  appendices  at  the  end  of  this  document.  A  block  flow  diagram  of  program 
usage  is  provided  in  Figure  3.1. 
1.      Program  File  thesis. h 

The  program  header  file  thesis. h  (Appendix  A)  contains  all  requisite  hbrary 
include  files,  C  manifest  constants,  and  MACRO  definitions  as  necessitated  by  main 
programs.  For  code  portability,  it  is  imperative  that  a  centralized  copy  of  all  common 
declarations  and  definitions  be  maintained  as  the  programs  evolve.  Therefore,  any 
modifications  to  common  parameters  can  be  performed  by  substitutions  within  this 
file. 

With  the  advent  of  ITEX/FCplus  software,  additional  include  files  were 
necessary  beyond  the  scope  of  the  original  header  file  listed  in  Reference  4.  For 
example,  ITEX/PC  libraries  require  the  use  of  the  include  file  itexpc.h  [Ref.  12:p. 
1-15].  This  requirement  has  been  superceded  by  the  ITEX/PCplus  software  which 
dictates  the  use  of  itexpfg.h  and  stdtyp.h  include  files.  Not  only  must  these  include 
files  be  listed  in  the  header  file,  but  to  be  utilized,  they  must  also  be  resident  in  the 
actual  program  directory  [Ref.  6:p.  1-9]. 

The  added  capability  of  dual  memory  required  that  many  of  the  calls  to  the 
software  libraries  be  altered.  Additional  manifest  constants  were  needed  for  memory 
specification.  They  included  MEMORY,  defining  the  type  of  frame  memory  con- 
figuration, and  constants  necessary  to  define  the  size  of  the  selected  frame  memory 

17 


speckle. c 


vir_8rry.c 


genfunc.c 


geofll.c 

Istol.c 

2sigm8.c 


threshll.c 


feel-id. c 


sizeit.c 


STATGRAPHICS 
SOFTWARE 


FILTERING 
ROUTINES 


IMAGE 
THRESHOLD 


FEATURE 
IDENTIFICATION 


FEATURE 
SIZING 


HISTOGRAM 
PRODUCTION 


Figure  3.1:  Flow  Diagram  of  Image  Processing  Programs. 
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configuration:  XSIZE,  YSIZE,  and  DEPTH.  XSIZE  and  YSIZE  indicate  respectively 
the  horizontal  and  vertical  coordinates  of  the  memory  size,  and  DEPTH  refers  to  the 
number  of  bits  assigned  per  pixel  in  frame  memory  [Ref.  6:p.  2-6]. 
2.      Program  File  geofiic 

Appendix  B  contains  the  listing  for  geofiic,  a  program  ba^ed  upon  a  geo- 
metric hulling  algorithm  introduced  in  References  9  and  10.  The  algorithm  is  applied 
to  the  image  and  its  complement  by  performing  comparisons  between  a  central  pixel 
value  and  adjacent  pixels  in  horizontal,  vertical  and  diagonal  directions.  The  original 
geometric  filter  routines  were  adapted  by  Edwards  [Ref.  1]  from  information  set  forth 
in  References  9  and  10;  and  Kaeser  successfully  adapted  these  to  the  C  language  [Ref. 
4].  Edwards'  FORTRAN  filtering  routine  was  capable  of  handling  only  a  quarter  of 
the  actual  image  due  to  memory  restrictions,  whereas  the  C  language  version  fully 
supported  the  entire  image.  A  large  part  of  the  C  version's  success  was  due  to  the 
utihzation  of  virtual  memory  arrays  for  data  processing.  These  arrays  are  stored  on 
disk,  but  they  are  accessed  as  though  they  are  stored  in  memory.  File  constraints 
within  the  associated  operating  system  po^e  the  only  limitation  upon  the  virtual  array 
size  [Ref.  13]. 

Few  modifications  were  performed  on  the  geofiic  program  to  make  it  com- 
patible with  the  ITEX/PCp/u5  software.  The  majority  of  the  changes  for  compatibil- 
ity encompassed  general  support  functions.  However,  much  of  the  C  language  code 
was  altered  to  attain  optimum  modularity.  This  filter  program  had  incorporated  the 
use  of  GO  TO  statements  to  facilitate  branching  between  the  filter  stepping  functions. 
Alterations  included  the  removal  of  these  subject  statements  and  the  formation  of  two 
recursive  functions.  These  recursive  functions  support  program  modularity  and  per- 
mit easier  modifications  for  future  updates. 
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One  common  modification  performed  on  geofil.c  and  its  other  main  pro- 
gram counterparts  involved  the  declaration  of  the  image  filename  and  comment  line 
(comhne).  The  header  file  used  by  ITEX/PC,  itexpc.h,  required  special  file  type  struc- 
tures used  in  conjunction  with  parameter  passing  for  SAVIM  and  READIM  functions. 
These  structures  consisted  of  strings  with  length  20  or  200  characters,  and  were  clas- 
sified as  typedef  LS20  or  LS200,  respectively  [Ref.  12:p.  1-19].  ITEX/PCp/u5  limited 
filename  to  20  characters  and  comment  lines  to  200  characters,  but  removed  the 
special  structure  requirement.  SAVIM  and  READIM  routines  now  require  that  the 
image  filename  and  comment  fine  be  classified  as  pointers  of  type  character  [Ref.  6]. 
Consequently,  all  program  files  were  altered  to  reflect  this  requirement. 
3.      Program  File  Istat.c 

The  program  Istat.c  (Appendix  C)  uses  local  statistical  methods  to  calcu- 
late a  weight,  k,  which  determines  the  gray  level  of  a  pixel  [Refs.  9  and  10].  The 
mean  and  variance  local  estimates  within  a  5x5  window  encompassing  the  selected 
pixel  determine  the  value  of  k.  The  k  value  then  specifies  the  new  gray  level  of  the 
pixel  in  question  [Ref.  l:p.  24].  Virtual  arrays  again  assist  in  the  data  manipulation 
and  storage  processes. 

Kaeser's  version  of  this  filtering  algorithm  required  that  the  user  input  the 
standard  deviation  of  the  image  directly  into  the  program  [Ref.  4].  The  operator  had 
to  obtain  this  value  from  the  histogram  function  of  the  ImageAction  software  [Ref. 
14:p,  11-9-3]  before  running  the  filter  program.  ImageActionp/u5  software  no  longer 
provides  the  image  mean  and  standard  deviation  calculations  [Ref.  5].  Consequently, 
a  new  function,  dev(stddev),  was  devised  to  calculate  the  image  standard  deviation 
and  input  it  directly  into  the  filter  routine  without  user  intervention.  This  new 
function  presents  a  vast  improvement  in  terms  of  user  convenience. 
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4.  Program  File  2sigma.c 

The  2sigma  filter  algorithm  provides  the  basis  for  2sigma.c  (Appendix  D) 
and  was  adapted  from  the  program  introduced  by  Edwards  [Ref.  1].  From  the 
standard  deviation  of  a  Gaussian  distribution,  it  is  known  that  95.5%  of  a  selected 
region's  pixels  will  be  within  two  standard  deviations  of  the  local  mean.  An  averaging 
formula  has  been  devised  in  which  pixels  meeting  the  "2sigma"  criteria  are  averaged 
for  a  given  pixel,  and  those  outside  the  range  are  rejected  [Ref.  10].  These  calculations 
are  performed  for  5x5  arrays  of  pixels.  The  program  2sigma.c  utilizes  virtual  arrays 
and  the  standard  deviation  function  as  well.  It  also  has  reduced  size  due  to  the 
removal  of  all  general  supporting  functions  into  the  program  file  genfunc.c. 

5.  Program  File  threshit.c 

Threshit.c  (Appendix  E)  produces  a  binary  version  of  the  filtered  image. 
The  program  operator  enters  a  threshold  value  which  sets  all  pixel  values  above  it  to 
WHITE  (level  255),  and  all  those  below  to  BLACK  (level  0).  The  user  then  views 
the  thresholded  image  and,  if  the  result  is  not  satisfactory,  is  given  the  option  to 
enter  a  new  threshold  limit.  Threshit.c  utilizes  a  function  named  threshit()  which  is 
physically  located  within  the  genfunc.c  program  file.  This  organization  was  due  to 
the  fact  that  other  program  files,  such  as  feat.id.c,  call  the  threshitf)  function. 

The  key  to  the  threshold  process  is  a  function  call  in  the  ITEX/PCplus  soft- 
ware package  called  THRESHOLD.  The  ITEX/PC  function  is  also  named  THRESH- 
OLD and  has  the  following  format: 

void  threshold (/owcw/,  highcut); 

The  lowcut  and  highcut  integer  values  defined  the  lower  and  upper  boundaries  of  the 
thresholded  region.  This  threshold  procedure  could  only  be  used  on  the  currently 
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selected  LUT  [Ref.  12:p.  7-27].  The  function  in  ITEX/PCp/u5  software  used  to  set  a 
selected  LUT  to  a  prescribed  threshold  is: 

void  threshold (^ro«p,  bank,  highcut,  lowcut); 

where, 

group  specifies  the  LUT:  RED,  GREEN,  BLUE,  or  INPUT, 
bank  specifies  the  output  LUT  used,  in  the  range  of  0  to  7, 
highcut  specifies  the  upper  bound  for  the  threshold,  and 
lowcut  specifies  the  lower  bound  for  the  threshold. 

The  threshit.c  program  utilizes  all  four  LUT  groups  and  the  0th  256-byte  bank  to 
specify  a  linear  transform  output.  This  specification  differs  significantly  from  the 
original  ITEX/PC  function  in  the  amount  of  versatility  granted  the  programmer. 
Any  of  the  LUTs  can  be  modified  with  this  threshold  command,  not  just  the  current 
LUT  [Ref.  6:p.  7-35].  The  lowcut  value  must  never  be  specified  as  higher  than  highcut 
value,  but  these  values  may  be  made  equivalent  to  display  one  intensity  at  level  255 
(WHITE). 

Once  the  user  is  satisfied  with  a  particular  image  threshold,  the  image  can 
be  mapped  to  memory.  This  procedure  involves  the  use  of  another  ITEX/PCp/us 
function: 

int  maplut (^rowp,  bank,  x,  y,  dx,  dy)\ 

M.A.PLUT  modifies  actual  memory  values  rather  than  just  LUT  values  as  THRESH- 
OLD does.  The  group  and  bank  variables  are  identical  to  those  used  in  the  THRESH- 
OLD function,  but  now  only  the  INPUT  LUT  is  modified  for  mapping  purposes  in 
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threshit.c.  The  quantities  x,  y,  dx,  and  dy  refer  to  the  horizontal  and  vertical  coor- 
dinates of  the  512x512  image  [Ref.  6:p.  7-17].  As  in  the  case  of  the  THRESHOLD 
function,  ITEX/PC  permits  the  user  to  transform  only  that  area  specified  by  the 
current  LUT.  The  variables  group  and  hank  are  not  defined  in  ITEX/PC  MAPLUT 
function  [Ref.  12: p.  7-14]. 

6.  Program  File  feat-id. c 

It  is  mandatory  that  an  image  is  thresholded  prior  to  initiating  the  program 
featJd.c  (Appendix  F).  For  this  reason,  the  user  is  given  the  option  within  featJd.c 
to  threshold  the  image  if  thresholding  was  not  completed  previously.  Conversion  to 
a  feature-labeled  image  requires  an  iterative  process,  whereby  the  image  is  evaluated 
row  by  row.  Once  this  feature  identification  process  is  complete,  the  image  can  be 
further  processed  for  specific  feature  data. 

The  featJd.c  program  necessitated  few  modifications  for  compatibihty  with 
ITEX/FCplus  software.  Again,  the  more  extensive  changes  were  performed  on  the 
general  support  functions.  FeatJd.c  used  the  ITEX/PCplus  functions  WPIXEL 
and  RPIXEL  extensively,  but  these  calls  retained  the  same  format  as  that  used  in 
ITEX/PC  software  [Refs.  6  and  12]. 

7.  Program  File  sizeit.c 

Program  sizeit.c  (Appendix  G)  creates  an  output  file  size.dat  which  con- 
tains the  dimensions  of  all  image  features  in  a  format  necessary  for  input  into  a  sta- 
tistical analysis  program.  Sizeit.c  relies  upon  pointer  structures  for  dynamic  memory 
allocation  of  processed  feature  data. 

A  conversion  factor,  SCALE-FACTOR,  figures  prominently  in  this  sizing 
algorithm.  Kaeser  determined  this  factor  to  be  3.7353  contingent  on  image  calibration 
[Ref.  4].  Upon  subsequent  runs  of  the  sizeit.c  program,  however,  a  SCALE_FACTOR 
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of  10.0000  was  necessary  to  achieve  compatible  results  with  Kaeser's  program  [Ref. 
4:p.  17].  The  constant  term  SCALE-FACTOR  is  defined  in  the  header  file  thesis.h 
and  can  be  easily  modified  if  required.  Bcised  upon  equipment  configuration  during 
hologram  processing,  the  subject  conversion  factor  is  still  questionable  in  terms  of 
accuracy.  For  this  reason,  an  additional  function  was  devised  to  determine  a  con- 
version factor  based  upon  user  input.  This  function  was  entitled  scale(),  and  its 
characteristics  are  discussed  in  the  section  relevant  to  genfunc.c. 

8.  Program  File  speckle. c 

To  calculate  the  value  of  speckle  index  (SI)  after  each  filter  iteration, 
speckle. c  (Appendix  H)  is  called  to  evaluate  filtering  effectiveness.  This  calculation 
is  performed  upon  specified  dimensions  of  the  filtered  image;  integer  values  for  im- 
age row  and  column  are  passed  to  this  function  from  the  filter  routines.  Values  for 
speckle  index  typically  range  from  a  maximum  value  of  1.0,  a  rare  occurrence,  to 
a  hypothetical  value  of  0  [Ref.  l:p.  20].  No  program  modifications  were  necessary 
for  ITEX/PCplus  compatibility.  Speckle  index  is  obtained  purely  from  statistical 
manipulations. 

9.  Program  File  vir.arry.c 

Vir_arry.c  (Appendix  I)  consists  of  a  group  of  functions  supporting  setup 
and  access  of  the  virtual  disk- based  array  stratagem  employed  primarily  by  the  filter 
routines.  The  virtual  array  functions  were  adapted  by  Kaeser  [Ref.  4]  from  Refer- 
ence 8.  These  arrays  reside  on  disk  and  are  accessed  as  though  they  are  located  in 
memory.  This  feature  overcomes  the  memory  limitation  normally  corresponding  to 
array  usage.  Through  the  use  of  MACROS  and  pointer  notation,  the  C  programming 
language  provides  the  means  of  managing  virtual  arrays  [Ref.  8].  Thesis.h  contains 
the  obligatory  MACROS  needed  for  structure  definition.  All  virtual  array  functions 
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mandated  the  use  of  long  integers  for  indexing  purposes  [Ref.  13].  Since  the  dimen- 
sions of  the  actual  image  are  defined  in  terms  of  integer  values,  a  conversion  interface 
problem  developed  in  all  the  filter  programs.  This  was  further  complicated  by  the 
fact  that  the  ITEX/FCplus  WPIXEL  function  references  datum  points  in  terms  of 
integer  values.  This  function  is  essential  for  "writing"  the  filtered  pixel  values  to 
memory  from  the  virtual  array  routines  [Ref.  6:p.  5-4].  With  the  use  of  the  C 
language  "(long)"  command,  the  integer  indices  referencing  pixel  values  within  the 
image  can  be  used  simultaneously  to  reference  values  within  the  virtual  arrays  with 
no  conversion  difficulty  [Ref.  15:p.  36]. 
10.      Program  File  genfunc.c 

The  program  file  genfunc.c  (Appendix  J)  contains  image  processing  general 
support  functions.  These  functions  maintain  routine  evolutions  occurring  within  the 
main  image  processing  programs.  The  following  sections  discuss  functions  that  were 
either  created  or  modified  due  to  the  advent  of  the  ITEX/PCp/u5  Hbraries. 

a.      Function  startit() 

This  function  makes  provisions  for  initial  Frame  Grabber  board  setup. 
All  the  main  function  files  with  the  exception  of  speckle,  c  and  vir.arry.c  access  this 
function.  Due  to  the  introduction  of  the  dual  memory  configuration,  new  ITEX/F Cplus 
calls  were  required  in  addition  to  what  Kaeser  utilized  [Ref.  4].  The  SETDIM  rou- 
tine determines  the  configuration  of  the  memory,  either  dual-store  or  single-store  [Ref. 
6:p.  2-6].  SELECT.MEM  determines  which  frame  memory  will  be  accessed  while  in 
the  dual-store  mode.  Frame  memory  choices  consist  of  MEM_A  or  MEM_B  as  de- 
fined in  the  itexpfg.h  header  file  [Ref.  6:p.  2-19].  DISPLAY _MEM  must  be  used  in 
conjunction  with  the  SELECT.MEM  routine.  It  specifies  which  memory  area  will  be 
displayed  and  gives  the  user  the  option  of  writing  to  one  image  area  while  displaying 
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the  other  [Ref.  6:p.  2-20].  For  the  current  apphcation,  MEM_A  is  both  the  accessed 
and  displayed  frame  memory. 

b.  Function  rtadit() 

The  image  is  read  from  disk  storage  and  placed  on  the  video  monitor 
via  the  readit()  function.  Any  comment  associated  with  the  image  when  it  was  stored 
will  be  displayed  at  that  time.  The  ITEX/PCp/u5  function  READIM  is  the  key  to 
this  procedure.  The  READIM  call  format  is  relatively  unchanged  from  the  ITEX/PC 
version.  If  an  error  occurs  when  a  file  is  being  read,  the  function  will  return  an  error 
code  to  identify  why  the  operation  failed.  The  itexpfg.h  header  file  defines  these  error 
code  values  [Ref.  6:p.  6-5].  The  standard  ITEX/PC  format  did  not  make  provisions 
for  these  values  in  its  associated  header  file,  itexpc.h;  therefore,  file  errors  are  referred 
to  numerically  [Ref.  12:p.  6-7].  Accordingly,  the  readit()  function  was  modified  to 
reflect  the  error  nomenclature. 

c.  Function  saveit() 

The  saveit()  function  allows  the  user  to  save  the  image  to  disk  memory 
via  the  ITEX/PCp/w5  function  SAVIM.  The  image  can  be  stored  using  compression; 
although  compression  may,  in  rare  cases,  cause  the  image  to  occupy  more  disk  space 
than  an  uncompressed  version.  If  compression  is  used,  the  file  will  be  automatically 
decompressed  when  read  from  the  disk.  As  in  the  case  of  the  READIM  function, 
SAVIM  returns  error  code  values  if  the  save  operation  is  unsuccessful.  Subsequent 
modifications  were  performed  to  the  saveit()  code  to  make  it  compatible  with  the 
new  documentation  [Refs  6  and  12]. 

d.  Function  dev(stddev) 

Both  program  files  Istat.c  and  2sigma.c  require  the  computation  of  the 
standard  deviation  of  the  image  for  their  respective  filtering  algorithms.  The  function 
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dev(stddev)  computes  this  value  through  an  iterative  process  utilizing  the  RPIXEL 
function.  Since  the  standard  deviation  is  computed  directly  from  the  image  without 
user  input,  there  is  a  slight  increase  in  filter  program  execution  times  associated  with 
the  implementation  of  this  function. 

e.      Function  scale  (factor) 

This  function  supports  the  sizeit.c  program  by  allowing  the  user  flex- 
ibility in  the  determination  of  a  scale  factor  used  in  the  sizing  algorithm.  The  user 
is  given  the  option  of  using  the  SCALE-FACTOR  as  delineated  in  the  header  file 
thesis. h^  or  of  entering  dimensions  of  a  calibration  object.  If  the  user  decides  upon 
the  latter,  the  length  of  the  object  must  be  entered  in  terms  of  pixels  and  inches. 
The  function  will  then  calculate  the  new  SCALE-FACTOR  and  pass  this  value  to  the 
sizeit.c  program. 

E.     SUMMARY 

The  modifications  and  improvements  performed  on  the  image  processing  pro- 
grams were  necessary  to  make  them  compatible  with  the  environment  established  by 
the  installment  of  the  PCVISIONp/us  Frame  Grabber  board  into  the  IBM  PC/AT. 
These  subject  programs  support  all  requirements  for  image  acquisition,  filtering, 
thresholding  and  particle  identification.  The  following  chapter  examines  the  per- 
formance and  verification  of  these  programs  at  length. 
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IV.  SOFTWARE  PERFORMANCE  ANALYSIS 

This  chapter  will  discuss  the  pertinent  details  concerning  software  performance. 
The  subjects  of  executable  program  size,  timing  analysis,  and  program  verification 
will  be  examined  at  length  in  the  following  sections.  Any  departures  from  results 
obtained  from  the  original  C  language  programs  will  be  appraised. 

A.     EXECUTABLE  PROGRAM  SIZE 

Table  4.1  shows  the  comparison  between  executable  program  size  of  ITEX/PC 
and  ITEX/PCplus  supported  software.  The  size  information  is  presented  in  terms  of 
kilobytes  (KB)  and  was  acquired  via  the  use  of  the  DOS  command,  DIR  [Ref.  4]. 
The  ITEX/PC- based  programs  written  by  Kaeser  [Ref.  4]  have  significantly  smaller 
executable  file  sizes.  This  is  attributed  to  the  fact  that  during  the  LINK  process, 
Kaeser  linked  main  programs  to  required  library  files  only.  It  was  not  necessary  to 
link  to  the  object  modules  of  the  genfunc.c^  speckle. c,  or  the  vir.arry.c  files  since  all 
necessary  functions  were  appended  to  the  end  of  the  main  program  files  [Ref.  4].  In 
addition,  the  ITEX/PC  function  library  differs  from  the  ITEX/PCp/w5  version,  and 
these  differences  will  also  impact  executable  file  size. 

TABLE  4.1:  Comparisons  of  Executable  Program  Size  (KB) 


FILENAME 

ITEX/PC 

ITEX/PCp/u5 

THRESHIT.EXE 

31,467 

63,101 

FEATJD.EXE 

32,509 

64,407 

SIZEIT.EXE 

46,281 

65,429 

2SIGMA.EXE 

49,483 

65,491 

LSTAT.EXE 

49,475 

65,499 

GEOFIL.EXE 

49,969 

66,317 
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TABLE  4.2:  Comparisons  of  Actual  Program  Size  (KB) 


FILENAME 

ITEX/PC 

ITEX/PCplus 

THRESHIT.C 

5,807 

1,336 

FEATJD.C 

10,130 

6,003 

SIZEIT.C 

11,011 

7,426 

2SIGMA.C 

14,623 

5,396 

LSTAT.C 

14,434 

5,396 

GEOFIL.C 

14,888 

7,916 

One  stated  objective  of  this  thesis  was  to  optimize  the  existing  C  language 
code.  Although  executable  program  size  comparisons  contradict  this  goal,  it  is  shown 
in  Table  4.2  that  actual  program  sizes  for  the  ITEX/PCplus  compatible  programs  are 
substantially  reduced  from  the  original  versions.  Two  linker  options,  /F  and  /PAC, 
were  initially  employed  to  achieve  more  compact  executable  files.  The  use  of  these 
options  should  produce  more  rapid  execution  of  files  and  shorter  program  load  times 
[Ref.  16:p.  185].  The  /E  option  had  been  utilized  by  Kaeser  in  Reference  4,  but 
this  file  packing  option  may  not  always  give  satisfactory  results.  It  can  occasionally 
increase  file  size,  and  this  obviously  reduces  savings  in  disk  space  [Ref.  17:p.  268]. 

Upon  Hnking  program  files  with  the  /PAC  and  /F  options,  it  was  found  that 
these  options  did  not  pose  an  improvement  over  the  /E  option.  In  fact,  executable 
file  sizes  were  increased.  For  example,  the  geofil.exe  file  was  73,132  KB  when  these 
options  were  used  in  the  LINK  process.  The  /E  option  produced  an  executable  file 
size  of  66,317  KB.  All  other  executable  files  were  affected  similarly.  Accordingly,  the 
/E  option  was  utilized  for  the  LINK  procedure.  Appendix  K  contains  the  listings  of 
the  batch  files  used  to  compile  and  fink  the  program  files. 
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TABLE  4.3:  Program  Execution  Times 


PROGRAM 

TIME(SEC) 

FEATURE  IDENTIFICATION 

809.82 

FEATURE  SIZING 

1780.03 

2SIGMA  FILTER 

760.55 

LOCAL  STATISTICS  FILTER 

1324.47 

GEOMETRIC  FILTER 

3403.13 

B.     PROGRAM  EXECUTION  TIMES 

The  jl7res4-img  library  image  was  used  to  test  program  execution  times.  This 
was  the  same  image  as  that  used  by  Kaeser  in  Reference  4,  depicting  an  Air  Force  Res- 
olution Target  (AFRT)  obscured  by  speckle  noise.  The  unfiltered  ima.ge  jl7res 4- in g 
is  depicted  in  Figure  4.1.  Testing  was  accomplished  by  the  use  of  a  timer  program, 
timer. c,  in  conjunction  with  a  DOS  input  file  to  eliminate  user  intervention.  Previous 
thesis  work  involved  time  testing  of  one  quarter  of  the  512x512  image.  MS  FORT 
versions  of  the  programs  could  effectively  handle  only  one  quadrant  of  the  full  screen 
image.  Consequently,  Kaeser  conducted  execution  time  comparisons  on  only  one 
quarter  of  the  image  [Ref.  4].  Due  to  the  banded  nature  of  the  ITEX/FCplus  frame 
memory,  quadrant  timing  tests  of  the  image  would  not  produce  results  that  could  be 
compared  with  the  original  findings  [Ref.  5];  therefore,  timing  analysis  of  the  pro- 
grams was  performed  on  the  entire  image.  The  resultant  execution  times  are  listed 
in  Table  4.3. 


C.     PROGRAM  VERIFICATION 

To  verify  the  results  obtained  from  the  ITEX/PCp/us-based  programs,  two  im- 
ages were  used  from  the  existing  image  library.  The  jllres^-img  was  again  used,  this 
time  to  determine  the  accuracy  of  the  filter  programs.  Both  feat-id. c  and  sizeit.c  were 
operationally  tested  with  the  image  lOxwgrid.img  derived  from  a  standard  calibration 
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Figure  4.1:  Unfiltered  jl7res4.img. 

image  [Ref.  2].  Resultant  data  acquired  from  these  images  was  compared  to  data  in 
Reference  4. 

Speckle  index  value  comparisons  were  used  to  determine  the  effectiveness  of  the 
filter  programs.  In  each  case,  the  respective  filter  was  run  through  six  iterations  of 
the  full  screen  image,  jl7res4-img.  Table  4.4  contains  the  tabulation  of  the  results 
acquired,  and  the  results  from  the  original  C  programs  are  listed  in  Table  4.5.  A  file 
comparison  was  performed  to  ensure  that  the  image  used  for  testing  was  identical 
to  the  full  screen  image  used  previously.  Comparison  of  the  data  obtained  from  the 
current  iterations  to  data  from  Kaeser's  programs  indicates  slight  variations  in  the 
speckle  index  values.  The  initial  speckle  index  values  for  the  two  cases  differ.  One 
possible  factor  attributing  to  this  difference  is  the  inconsistency  of  integer  conversions 
between  the  C  language  programs  [Ref.  4].  The  more  significant  measure  of  filter 
performance  is  the  relative  decrease  in  speckle  index  upon  subsequent  filter  iterations. 
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TABLE  4.4:    Speckle  Index  Reduction  Results  Obtained  from  Modified 
Programs 


NUM 

2  SIGMA 

L  STAT. 

GEOMETRIC 

0 

0.293046 

0.293046 

0.293046 

1 

0.196282 

0.180061 

0.218982 

2 

0.144286 

0.125157 

0.174153 

3 

0.115754 

0.099063 

0.143760 

4 

0.098420 

0.084035 

0.121202 

5 

0.087068 

0.074063 

0.103756 

6 

0.079042 

0.066919 

0.089894 

TABLE  4.5:  Speckle  Index  Reduction  Results  (Ref.  4) 


NUM 

2  SIGMA 

L  STAT. 

GEOMETRIC 

0 

0.304005 

0.304005 

0.304005 

1 

0.153767 

0.194121 

0.224875 

2 

0.098967 

0.151951 

0.178222 

3 

0.074415 

0.129769 

0.146152 

4 

0.060615 

0.115938 

0.122347 

5 

0.051691 

0.106372 

0.103983 

6 

0.045354 

0.099367 

0.089425 

Therefore,  the  dissimilarities  of  initial  speckle  index  value  are  not  detrimental  to  filter 
function,  and  all  filter  programs  produce  desirable  reduction  of  speckle.  Consequently, 
all  evidence  indicates  that  the  ITEX/PCp/us-based  filter  programs  function  correctly. 
Figure  4.2  illustrates  jllres^.img  after  six  iterations  of  the  geometric  filter. 

The  results  obtained  from  the  feat  Ad. c  and  sizeit.c  programs  were  identical  to 
those  obtained  from  Kaeser's  programs  [Ref.  4].  The  SCALE-FACTOR  entry  in  the 
thesis. h  header  file  had  to  be  modified  in  order  for  the  output  files  from  the  respective 
sizing  programs  to  be  congruent.  Otherwise,  the  program  files  produced  the  same 
outputs  as  the  original  C  language  programs. 
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Figure  4.2:  J17res4.img  Image  After  Six  Iterations  of  the  Geometric  Filt 
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V.  CONCLUSIONS 

A.  SUMMARY  OF  FINDINGS 

This  thesis  accomplished  all  major  objectives  as  set  forth  in  the  introductory 
chapter.  No  major  problems  were  anticipated  with  the  integration  of  the  ITEX/PCp/M5 
libraries  and  the  COMPAQ  386/20  computer  system.  Due  to  a  memory  conflict  cre- 
ated within  the  COMPAQ,  the  PCVISI0Np/u5  had  to  be  installed  into  the  IBM 
PC/AT  to  complete  C  program  modifications.  Therefore,  the  bulk  of  this  research 
was  conducted  on  the  PC/AT  used  in  previous  work.  Although  several  of  the  pro- 
grams produced  in  Kaeser's  research  [Ref.  4]  were  modified  with  no  major  difficulties, 
type  definitions  created  major  faults  in  the  filter  routines  which  had  to  be  resolved. 
It  was  the  goal  of  this  thesis  to  improve  program  portability  and  optimization.  For 
the  most  part,  this  goal  was  achieved.  Executable  program  size  did  increase,  but  this 
was  the  result  of  program  structure  changes  and  library  modifications  introduced  by 
the  ITEX/Pcp/u5  system  upgrade. 

B.  FUTURE  WORK 

Further  verification  testing  is  necessary  to  suflficiently  demonstrate  the  operation 
of  the  C  language  programs.  A  minimal  amount  of  experimentation  was  accomphshed 
to  compare  the  acquired  results  with  data  obtained  from  the  previous  version  of  the 
programs.  Further  testing  is  mandatory  for  complete  program  validation  and  should 
be  pursued  in  subsequent  thesis  research  prior  to  porting  the  programs  over  to  other 
computer  systems.  Some  effort  was  expended  on  the  application  of  the  NDP  C-386 
compiler  to  eliminate  the  need  for  virtual  arrays  for  data  processing.  These  efforts 
were  largely  unsuccessful  due  to  the  incompatibility  of  the  YTEX/FCplus  libraries 
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with  the  subject  compiler  system.  This  area  of  research  could  be  further  investigated 
for  future  thesis  work.  The  following  topics  may  also  be  considered  for  future  work: 

•  Utilization  of  Macintosh  board  for  image  processing. 

•  Measurement  of  image  particles  utilizing  a  scanning  electron  microscope  for 
sizing  routine  validation. 

•  Further  validation  of  image  measurements  using  new  caHbration  objects. 

•  Integration  of  SUN  workstations  for  image  processing. 
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APPENDIX  A:  PROGRAM  HEADER  FILE: 

thesis. h 


/*      PURPOSE: 

To  completely  list  all   necesseury  include  files,  mzmifest  constemts, 
and   MACRO  definitions  required  by  all  the  progreims  sithin  this 
package.   Any  changes  to  constants  in  thesis. h  cill  guaxaintee 
redefinition  of  those  paxauneters  throughout  the  main  files. 

*/ 

/*     Include  files  for  use  with  ITEXPC  programs 
*/ 

#include  <stdio.h> 
#include  <stdlib.h> 
#include  <math.h> 
#include  <malloc.h> 
#include  <graph.h> 
#include  <conio.h> 
#include  "itexpfg.h" 
#include  "stdtyp.h" 


/*   Progreim  constants  used  in  the  C  program  modules  of  the  ITEXPC 

Optical  System  used  for  Rocket  motor  Hologram  amalysis.   */ 

/*  Intial  ITEXPC  Boaxd  Jumper  Settings   */ 

#define  MEMBASE    OxDOOOOL   /*  Itex  board  base  memory  start  address  */ 
#define  REGBASE    0x300      /*  Itex  boaxd  base  register  start  address  */ 
#define  MEMORY     DUAL       /*  Itex  board  memory  type  */ 

/*  ITEXPC  initial  AOI  (area  of  interest)  settings  */ 


#def ine 

IXS 

0 

#def ine 

lYS 

0 

#def ine 

NROW 

480 

#def ine 

NCOL 

512 

#def ine 

COLOR 

150 

#def ine 

XSIZE 

512 

#def ine 

YSIZE 

512 

#def ine 

DEPTH 

8 

/*  ItexPC  LUT  (Look  Up  Table)  Variables      */ 
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/*  Threshold  Limits  ♦/ 

#define  LOWEST    0  /*  Equates  to  Black  lor  lowcut  value  */ 

#define  HIGHEST   255        /*  Equates  to  White  for  highcut  value  */ 

#define  BLACK     0 

#define  WHITE     255 

/♦  Filter  array  sizes  and  limits  */ 

#define  SHUM  25 

#define  HUM  9 

#define  HIGH  254 

#define  LOW  0 

/*  Sizing  magnification  factor  used  in  sizing  program  */ 

#define  SCALE.FACTOR  10.000000   /*  Based  on  (  317.5  /  85  )   */ 

/*  Virtual  airray  Header  File  Definitions   */ 
#define  header        7 

/*  Virtual  Array  Control  Block  typedef  */ 

typedef  struct  { 

FILE  *file;  /*  pointer  to  file  control  block  */ 

long  size;  /*  number  or  array  elements  in  file  */ 

int  elsize;  /♦  number  or  bytes  in  each  element  */ 

char  *buffer;  /*  pointer  to  array  buffer  */ 

int  buf_elsi2e;  /*  size  of  element  in  buffer  including  index  */ 

int  buf_size;  /*  number  of  elements  in  buffer  */ 

char  *blcink_rec;  /*  pointer  to  initialization  record  */ 

/*  used  for  extending  file  */ 

} 

VACB  ;  /*  Virtual  Array  control  block  type  najne  */ 

/*  Virtual  Array  Access  Prototypes  ♦/ 

int  init_v_axray(cheLr  *f  ilename.int  rec_size,chaT  f  ilchar)  ; 

VACB  *open_v_array(chELr  *f  ilename,  int  buff  er_size)  ; 

void  close_v_array(VACB  *v_array) ; 

void  *access_v_rec(VACB  *v_array .long  index); 

/*  Virtual  Array  Access  Macros  */ 

#define  VREC(y)  ((items  *)access_v_rec(item_eLrray  ,y) ) 
#define  gdesc(y)  VREC(y)->v_gdesc 
#define  G(x,y)  VREC(y)->v_gdesc[x] 

/♦  Virtual  Array  element  structxire  typedef  */ 
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typedef  struct  { 

unsigned  char  v_gdesc[NCOL] ; 
} 
items; 
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APPENDIX  B:  PROGRAM  FILE:  geofil.c 


/♦PURPOSE   :  Provides  for  filtering  an  imaige  using  the  geometric  filter 
algorithm.   This  program  processes  the  image  and  provides 
for  the  operator  to  select  the  n\unber  of  iterations. 

The  Eorea  of  filtering  for  the  image  is  controlled  by  the 
operator  by  keyboaird  input  of  number  of  ROWS  amd  COLS  to 
process.   Progreun  will  sdlov  for  a  Maximum  input  of  480  rows 
and  512  cols  with  a  minimum  of  1  row  sind  1  column. 

The  resulting  filtered  image  pixel  value  is  stored  on 
disk  in  a  virtual  array,  that  is  written  to  the  image 
processing  screen  when  eill  ceilculations  axe   complete. 

Totcil  calculation  time  for  a  512x480  image  is  3403.13  sec. 

DISK  SPACE  REQUIRED  for  virtual  array  storage  is  250  KB  minimum 
available  on  the  disk  at  staxt  of  the  routine. 

*/ 
#include  "thesis. h" 

char         filename [20] ,    comline[200] ; 
int  srow,    scol,    times,    row,    col; 

mainO 


int  flag; 

extern  int   row,  col,  srow,  scol,  times; 

printf("\n\tReady  to  Load  IMAGE  FILE?... Yes  (y)  /  No  (n)  "); 

flag=getch(); 

if (flag  ==  'Y'l I  flag  ==  'y') 

{ 

startitO  ; 

readitO; 
} 
if(  flag  ==   'N'  II  flag  ==   'n'  )  return(O) ; 

geti("\n\nEnter  niimber  of  rows  (Max=480)  to  use  in  SPECKLE  CALC  :  ".ftsrow); 
geti("\nEnter  number  of  columns  (Max=512)  to  us  in  SPECKLE  CALC  :   ",&scol); 
geti("\n\n\nEnter  number  of  filter  iterations  to  do. . . ".fttimes) ; 
geti("\nEnter  the  number  of  IMAGE  rows  (MAX  =  480)  to  FILTER  :  ",&row); 
geti("\nEnter  the  number  of  IMAGE  cols  (MAX  =  512)  to  FILTER  :  ",&col); 
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sBitch(times) 
{ 

case  0: 

speckleO ; 

printf ("\nSave  image  to  Disk  File?.. Yes  (y)  /  Ho  (n)  "); 

llag=getch(); 

iKflag  ==    'Y'll   Hag  ==    'y')   saveitO; 

exit(l): 

breaJc; 

case      1: 

break; 


> 

do 


default : 

breaJc; 


speckleO  ; 
geofilO; 


while(times   ==  0) ; 

printf ("\nSave   image  to  Disk  File?.. Yes    (y)   /  No   (n)    "); 

flag=getch() ; 

if  (flag  ==  'Y'll  flag  ==  'y')  saveitO; 


/*  geometric  filter  algorithm  */ 


geof il() 


VACB  ♦item_axray ; 

extern  int  row,  col,  times; 

int  fl,  f2,  gl,  g2,  g3,  X,  y; 

int  a,  b,  c,  d,  pixel,  calls; 

/♦  create  a  virtual  axray  for  the  array  of  filtered 
pixel  values  the  size  of  the  image  (  480x512  )    */ 

init_v_array (" ITEMS. VAR",sizeof (items) .NULL) ; 

/*  open  the  virtual  axray,  reserve  buffer  space  for  100  elements  */ 

item_array  =  open_v_array("ITEHS.VAR" , 100) ; 

/*  stairt  local  statistic  filter  routine   */ 

if(  row  >  NROW  )  row  =  NROW; 
if(  col  >  NCOL  )  col  =  NCOL; 
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lor(  y  =  1  ;  y  <  row  ;  y++  )      /*  Zero  the  initial  array  */ 
{ 

for(  X  =  1  ;  X  <  col  ;  x++  ) 
{ 

G(((long)x).((long)y))  =  HULL; 
} 
} 

for(calls  =  times;  calls  !=  0;  calls — ) 
{ 

c  =  4; 

printf  ("\n\n\tFilter  Running  — >  '/.d  Runs  after  this  one !  . .  .\n",calls-l)  ; 

lor(  ;  c  ;  ) 

{ 

8witch(c) 

{ 

case  4  : 

a=l;  b=0;  c=3;  d=l; 
printf ("\nc  =  3  "); 
tl(a,b,c,d,itein_array) ; 
brecik; 

case   3: 

a=0; 

b=li 

c=2; 

printf ("\nc  =  3  "); 

tl(a,b,c,d,item_array) ; 

breeA; 

case   2: 

a=li 

b=li 

c=li 

printf ("\nc  =  2  "); 

tl(a,b,c,d,item_2Lrray) ; 

breaik; 

case    1: 

a=l; 

b=-l; 

c=0; 

printf ("\nc  =  1  "); 

tl(a,b,c,d,item_aLrray) ; 

printf ("\nc  =  0  END"); 

speckleO  ; 

breaik; 
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default    :  breeik; 
} 
} 
} 

/*  close  the  virtual  array  */ 
close_v_aLrray( item. array) ; 
} 


/*  subroutine  lor  geolil  lor  steps  1  and  2  */ 

tl(a,b,c,d,itein_jurray) 
int    a,  b,  c,  d; 
VACB   *item_eu:Tay; 
{ 

int  11,  12,  gl.  g2,  g3; 

int         pixel,  x,  y; 
long         xl,  yl; 
extern  int   row, col, times; 

lor(   y  =   1    ;    y  <  row-l    ;    y++   ) 
{ 

lor(   X  =   1    ;    X  <  col-1    ;    x++   ) 
{ 

xl   =    ((long)  x) ; 

yl   =   ((long)  y); 

11  =  rpixeK  x,y   ) ; 

12  =  rpixeK  x-a,y-b   ); 

G(xl,yl)=  max(   11,   min(  12-1,    11+1    )    ); 
} 
} 

printK"      step   1"); 
lor(   y  =   1    ;    y  <  row-l    ;    y++    ) 
{ 

lor(   X  =   1    ;    X  <  col-1    ;    x++   ) 
{ 

xl  =  ((long)  x); 
yl  =  ((long)  y); 
gl  =     G(  xl.yl); 

g2  =   G( ((long) (x-a)),( (long) (y-b))); 
g3  =   G(((long)(x+a)).((long)(y+b))); 
pixel  =  max(  gl,  min(  g2,  min(  gl+1,  g3+l  )  )  ) 
wpixel(x,y, pixel) ; 
} 
} 

printK"   step  2") ; 
iK  d==l  ) 
{ 

a=  -a; 
b=  -b; 
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d=  0; 

tl(a,b,c,d,item_array) ; 
> 

else  if(  d==0  ) 
{ 

d=l; 

t2(a,b,c,d,item_axray) ; 
return ; 
} 


/♦  subroutine  lor  geolil  lor  steps  3  aaid  4  */ 
t2(a,b,c,d,item_eLrray) 
int  a,  b,  c,  d; 

VACB         *item_eirray; 


int  fl,  12,  gl,  g2,  g3; 

int  X,  y,  pixel; 

long  xl,  yl; 

extern  int  row,  col,  times; 

lor(  y  =  1  ;  y  <  row-l  ;  y++  ) 
{ 

lor(  X  =  1  ;  X  <  col-1  ;  x++  ) 
{ 

xl  =  ((long)  x); 

yl  =  ((long)  y); 

11  =  rpixel(  x,y  ) ; 

12  =  rpixel(  x-a,y-b  ); 

G(xl,yl)=  min(  11,  max(  12-1,  11+1  )  ); 
} 
} 

printK"   step  3") ; 
lor(  y  =  1  ;  y  <  row-l  ;  y++  ) 
{ 

lor(  X  =  1  ;  x  <  col-1  ;  x++  ) 
{ 

xl  =  ((long)  x); 
yl  =  ((long)  y); 
gl  =  G(  xl.yl); 

g2  =     G(((long)(x-a)),((long)(y-b))); 
g3  =     G(((long)(x+a)),((long)(y+b))); 
pixel  =  inax(  gl,   min(   g2,    min(   gl  +  1,   g3+l    )    )    ); 
wpixeKx.y, pixel); 
} 
} 

printl("      step  4") ; 
il(  d==l    ) 
{ 
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a=  -a; 
b=  -b; 
d=  0; 
t2(a,b,c,d,item_aurray) ; 


else   il(   d==0   ) 
{ 

d  =  1; 

return; 
> 
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APPENDIX  C:  PROGRAM  FILE:  Istat.c 


/♦PURPOSE  :  Provides  lor  filtering  an   imeige  using  the  local  statistical 

algorithm.  This  progrein  processes  the  image  by  using  a  local 
5x5  airray  of  pixels  to  calculate  a  statistical  value  for  the 
central  pixel  of  the  local  airray. 

The  eirea  of  filtering  for  the  image  is  controlled  by  the 
operator  by  keyboeird  input  of  number  of  ROWS  and  COLS  to 
process .  Program  will  eillov  for  a  Neucimum  input  of  480  rows 
and  512  cols  with  a  minimum  of  1  row  and  1  column.   Program 
also  requires  the  image  Std.  Dev.  as  calculated  by  the  genf unc. c 
Program.   The  resulting  filtered  image  pixel  value  is  stored  on 
disk  in  a  virtual  zorray,  that  is  written  to  the  image  processing 
screen  when  all  calculations  axe   complete. 

Total  calculation  time  for  a  512x480  image  is  1324.47  sec. 

DISK  SPACE  REQUIRED  for  virtual  array  storage  is  250  FB 
minimum  available  on  the  disk  at  the  start  of  the  routine. 
*/ 

#include  "thesis. h" 

char    filename [20] ,  comline[200] ; 
int     row,  col,  srow,  scol,  times; 
float   devl; 

mainO 

{ 

int  flag; 

extern  int    row,  col,  times; 

printf ("\n\tReady  to  Load  IMAGE  FILE?... Yes  (y)  /  No  (n)  ") ; 

flag=getch(); 

if (flag  ==  'Y'l I  flag  ==  'y') 

{ 

stctrtitO  ; 

readitO  ; 
} 
if (flag  ==  'N'  II  flag  ==  'n')  return(O); 

geti("\n\nEnter  number  of  rows  (Max=480)  to  use  in  SPECKLE  CALC  :  ".ftsrow); 
geti("\nEnter  number  of  columns  (Max=512)  to  us  in  SPECKLE  CALC  :   ",&scol); 
geti("\n\n\nEnter  number  of  filter  iterations  to  do. . . " ,&times) ; 
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geti("\nEnter  the  number  ol  IMAGE  rows  (MAX  =  480)  to  FILTER  :  ".krow); 
geti("\nEnter  the  number  ol  IMAGE  cols  (MAX  =  512)  to  FILTER  :  ".ftcol); 
dev(ftdevl) ; 
printl("\nThe  standard  deviation  lor  this  image  is  y,.2f\n"  ,devl) ; 


do 
■C 


speckle 0 ; 
IstatO; 

vhile(times   ==  0); 

printl("\nSave  image  to  Disk  File?.. Yes    (y)   /  lo   (n)   "); 

llag=getch(); 

iKflag  ==  'Y'll  Hag  ==  'y')  saveitO; 


/♦ 

IstatO 


.*/ 


VACB  *item_array; 

long  xl.yl; 

extern  float  devl; 

extern  int  row,  col,  times,  srow,  scol; 

register  int  m,  n,  nn; 

int  pixel, a, b,x,y .calls; 

float  suml,  sum2,  svar,  Ivar,  tveir; 

float  Imean,  lmean2,  stddev,  stddev2,  k; 

int  IdataCSNUM];   /*  SNUM  defined  in  Header  file  */ 

/*  create  a  virtual  eirray  for  the  eirray  of  filtered 
pixel  values  the  size  of  the  image  (  480x512  )    */ 

init_v_array("ITEMS.VAR",sizeof (items) .NULL) ; 

/*  open  the  virtual  array,  reserve  buffer  space  for  100  elements  */ 

item_array  =  open_v_array("ITEMS. VAR", 100) ; 

/*  stcirt  local  statistic  filter  routine   */ 

if(  row  >  NROW  )  row  =  NROW; 
if(  col  >  NCOL  )  col  =  NCOL; 


stddev  =  devl  /  255.0  ; 
stddev2  =  stddev  *  stddev  ; 


/*  Normedize  deviation  value  */ 


for(calls  =  times;  calls  !=  0;  calls — )   { 

printf  ("\n\n\tFilter  Running  — >  '/.d  Runs  after  this  one!  .  .  .\n"  ,calls-l) ; 
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lor  (  y  =  2  ;  y  <  row  -  2  ;  y++) 
{ 

lor  (  X  =  2  ;  X  <  col  -  2  ;  x++) 

{ 

pixel  =  rpixeK  x,y  ); 
suml  =  0; 
8uin2  =  0; 
nn  =  0; 

lor(  n=y-2;n<y+3;   n++   ) 

{ 

lor(  in  =  x-2    ;bi<x  +  3    ;   in++   ) 
{ 

ldata[nn]  =  rpiiel(  m,n  ); 

suml  +=  IdataCnn] ; 

nn++; 
} 


Imeaji  =     suml   /  SNUH      ; 

lme8Ln2  =   Imean  *  Imeaji   ; 

lor(   nn  =  0   ;    nn  <  SNUM    ;    nn++   ) 

{ 

svax  =   (ldata[nn]    -  ImeaJi   )    *    (   ldata[nn]    -  Imeem   ); 

sum2  +=  svar; 
} 

Ivar  =  sum2  /  SNUM  ; 

tvar  =  labs(  (  (Ivar  +  lmean2)  /  (stddev2  +  1)  )  -  lmean2  )  ; 
k  =  tvar  /  (  (  stddev2  *  lmeaji2  )  +  tvar  )  ; 
xl  =  ((long)  x) ; 
yl  =  ((long)  y) ; 

G(  xl.yl  )  =  (int)  (  Imeaji  +  k  *  (  pixel  -  Imeaji  )  )  ; 
} 
> 

/*  write  contents  of  the  G  array  to  image  screen  ♦/ 

print! ("\n\a\tWr it ing  filtered  image  to  screen..."); 

lor(  b  =  2  ;  b  <  row  -  2  ;  b++  ) 
{ 

for(  a  =  2  ;  a  <  col  -  2  ;  a++  ) 

{ 

wpixeK  a,b,G(((long)a),((long)b))); 

} 
} 

speckle  0 ; 
} 

/*  close  the  virtual  airray  */ 
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close_v_array(item_array) ; 
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APPENDIX  D:  PROGRAM  FILE:  2sigma.c 


/♦PURPOSE   :  Provides  for  filtering  an  imaige  using  a  2-sigma  statistical. 

eilgorithm.   This  program  processes  the  image  by  using  a  local 
5x5  array  of  pixels  to  calculate  a  statistical  value  for  the 
central  pixel  of  the  local  arr&y .     The  algorithm  vill  eliminate 
pixels  from  the  summation  that  axe   greater  than  2-Sigma  in 
value  from  that  of  the  centred  pixel. 

The  aurea  of  filtering  for  the  image  is  controlled  by  the 
operator  by  keyboaurd  input  of  number  of  ROWS  amd  COLS  to 
process.   Program  will  adlow  for  a  Maximum  input  of  480  rows 
and  512  cols  with  a  Minimum  of  1  row  and  1  coliunn.   Program 
also  requires  the  image  Stamdaird  deviation  as  obtained  using 
genfunc.c  prograun.   The  resulting  filtered  pixel  value  is 
stored  on  disk  in  a  virtual  aorray  amd  written  to  screen  later. 

Totail  calculation  time  for  a  512x480  image  is  760.55  sec 

DISK  SPACE  REQUIRED  for  virtual  array  storage  is  250  KB 
minimum  available  on  the  disk  at  stairt  of  the  routine. 

*/ 
#include  "thesis. h" 

char  filename [20] ,  comline[200] ; 

int  row,  col.t  imes,  srow,  scol; 

float  devl; 

mainO 


int  flag; 

extern  int     row,  col,  times; 

printf("\n\tReady  to  Load  IMAGE  FILE  from  disk?... Yes  (y)  /  No  (n)  "); 

flag=getch(); 

if (flag  ==  'Y'l I  flag  ==  'y') 

{ 

startitO  ; 

readitO  ; 
} 

if (flag  ==  'N'  II  flag  ==  'n')  return(O) ; 

geti("\n\nEnter  number  of  rows  (Max=480)  to  use  in  SPECKLE  CALC  :  ".ftsrow); 
geti("\nEnter  number  of  columns  (Max=512)  to  us  in  SPECKLE  CALC  :   ",&scol); 
geti("\n\n\nEnter  number  of  filter  iterations  to  do. . . " .fttimes) ; 
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geti("\nEnter  the  niunber  ot   IMAGE  rows    (MAX  =  480)   to  FILTER   :    ",»row); 

getiC'XnEnter  the  number  ol   IMAGE  cols    (MAX  =  512)  to  FILTER   :    ",»col); 

dev(*devl) ; 

printf ("\n\tThe  steoidard  deviation  for  this  image  is:     %.2f\n",devl) ; 

do 

{ 

speckleO ; 

sigmaO; 
} 
Bhile(  times  ==  0  ) ; 

printf ("\nSave  image  to  Disk  File?.. Yes  (y)  /  lo  (n)  ");    ^ 

flag=getch(); 

if  (flag  ==  'Y'll  flag  ==  'y')  saveitO; 


/^Hf**^:^i^;^iH^:^:^^iitHHHL:t:*t^:IHli;^^^;^t*************************************r)HL/ 


sigmaO 


VACB  *item_array; 

extern  int  row,  col,  times,  srow,  scol; 

extern  float  devl; 

long  xl,  yl; 

register  int  m,  n,  nn; 

int  pixel,  s\iml ,  sum2,  delta,  high; 

int  low,  a,  b,  x,  y,  calls; 

float  hvcir,  Ivar,  stddev; 

int  ldata[SNUM];         /*  SNUM  defined  in  Header  file  */ 

/♦  create  a  virtual  axray  for  the  array  of  filtered 

pixel  values  the  size  of  the  image  (  512x480  )        */ 

init_v_array("ITEMS.VAR",sizeof (items) ,NULL) ; 

/*  open  the  virtual  axray,  reserve  buffer  space  for  100  elements  */ 

item_array  =  open_v_array ( "ITEMS. VAR", 100) ; 

/*  staxt  local  statistic  filter  routine   */ 

if(  row  >  NROW  )  row  =  NROW; 
if(  col  >  NCOL  )  col  =  NCOL; 

stddev  =  devl  /  255.0  ;         /*  Normalize  deviation  value  */ 
hvar  =  1.  +2.  *  stddev  ;       /♦  find  high  and  low  sigma  */ 
Ivsir  =  1.  -  2.  *  stddev  ; 

for  (calls  =  times;  ceills  !=  0;  calls — )   { 

printf  ("\n\n\tFilter  running  — >  '/,d  Runs  after  this  one !  .  .  .\n",calls-l)  ; 
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lor  (  y  =  2   ;    y  <  row  -  2   ;   y++) 
{ 

for  (  X  =  2   ;    X  <  col  -  2    ;    x++) 

{ 

pixel   =  rpixelC   x,y   ); 

suml   =0; 

sum2  =  0; 

nn  =  0; 

delta  =  0; 

high  =   (int)(hvax  *  pixel); 

low  =   (int)(lvau:  *  pixel); 

lor(  n=y-2;n<y+3;   n++   ) 
{ 

for(  in  =  x-2;m<x  +   3;   ni++   ) 
{ 

xl   =    ((long)   x); 

yl  =   ((long)   y) ; 

IdataEnn]    =  rpixeK  m.n  ); 

if(    (   low  <=  IdataCnn]    )   ft&    (   ldata[nn]    <=  high   )    ) 

■C 

suml  +=  ldata[nn] ; 
delta++; 
> 
nn++; 


if(delta  <=  2)    /*  correct  shot  noise  —  4  neighbor  average  */ 
{ 

sum2  =  (  rpixel(x,y-l)+rpixel(x,y+l)+rpixel(x-l,y)+rpixel(x+l,y)  ); 

G(   xl.yl    )    =   (pixel   +  s\un2)    /   5; 

continue; 
> 

G(  xl.yl  )   =  sural  /  delta   ; 
} 
} 

printf ("\n\tWriting  filtered  image  to  screen...");  *'■ 

for(  b  =  2  ;  b  <  row  -  2  ;  b++  ) 
{ 

for(  a  =  2  ;  a  <  col  -  2  ;  a++  ) 

wpixeK  a,b,G(((long)a),((long)b))  ); 
} 

speckleO  ; 
} 

/*  close  the  virtual  array  */ 
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close_v_array(item_array) ; 
} 
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APPENDIX  E:  PROGRAM  FILE:  threshit.c 


/♦PURPOSE   :  This  program  thresholds  the  image  on  screen  by  talking 
an  operator  input  value  aind  forcing  all  image  pixel 
values  above  the  threshold  value  to  the  value  of  BLACK 
emd  all  those  below  the  value  to  the  WHITE  value. 

(  Background  ==  WHITE  ;  Feature  ==  BLACK  ) 


*/ 
#include  "thesis. h" 

char  filename  [20]  ,  coinline[200]  ; 

int  srow,  scol; 


mainO 
■C 


int  flag; 

printf ("\nThis  prograon  will  threshold  the  input  image  desired."); 

printf ("\n\n\tReady  to  Load  Image? .. .Yes  (y)  /  No  (n)  "); 

flag=getch(); 

if (flag  ==  'Y'l I  flag  ==  'y') 

{ 

stairtitO  ; 

readitO  ; 
} 

if (flag  ==  'N' I  I  flag  ==  'n') 
{ 

initializeO  ; 

return(O) ; 
}  ,, 

printf ("\n\n\tTHRESHOLD  the  image?... Yes  (y)  /  No  (n)  "); 

flag  =  getchO  ; 

if  (flag  ==  'YMI  flag  ==  'y')  threshitO; 
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APPENDIX  F:  PROGRAM  HEADER  FILE: 

feat  Jd.c 


/♦PURPOSE   :  Labels  and  identifies  each  feature  in  an   image  on  the  ITEX 
system.   Reads  pixel-by-pixel  and  groups  the  objects  by 
assigning  a  unique  ID  number  to  each  feature  or  object  so 
that  they  can   be  processed  by  other  routines. 

The  ID  routine  requires  a  thresholded  on  image  screen  image 
or  input  of  a  previously  saved  thresholded  image  from  disk  or 
imput  from  disk  and  then  thresholding.  This  progreun  module 
includes  an  optional  call  to  threshold  if  desired. 

(  Background  ==  WHITE  ;  Feature  ==  BLACK  ) 

*/ 
#include  "thesis. h" 

char    f  ileneime[20]  ,  comline[200]  ; 
int     srow,  scol; 


mainO 


int  flag; 

printf("\n\aYou  MUST  USE  A  THRESHOLDED  IMAGE  for  this  Program!!"); 

printf ("\n\n\tReady  to  Load  image  from  disk?... Yes  (y)  /  No  (n)  "); 

flag=getch(); 

if (flag  ==  'Y' I  I  flag  ==  'y') 

■c 

stairtitO  ; 

readitO  ; 
} 
if (flag  ==   'N'  II  flag  ==   'n'  )  return(O); 

printf ("\n\n\tNeed  to  THRESHOLD  the  image?... Yes  (y)  /  No  (n)  "); 

flag  =  getchO  ; 

if  (flag  ==  'Y'll  flag  ==  'y')  threshO; 

feat_id() ; 

printf ("\n\nSave  image  to  Disk  File?.. Yes  (y)  /  No  (n)  ") ; 

flag=getch(); 

if  (flag  ==  'Y'll  flag  ==  'y')  saveitO; 
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/*  image  feature  identification  algorithm  */ 

feat.idO 
■C 

register  int    x,  y; 

int  a,  xl,  yl,  x2,  y2,  x3,  y3,  x4,  y4; 

int  fid  =  0; 

int  gid  =  0; 

int  temp_fid  =  0; 

int  temp_gid  =  0; 

int  maxfl,  nl ,  n2,  nla,  nib,  nlc,  kk,  kkl,  nnl; 

printf ("\n\n\tSTEP  ONE  in  Progress..."); 

for(  y  =  1  ;  y  <  NROW  ;  y++  ) 

for(  X  =  0  ;  X  <  NCOL  ;  x++  ) 
{ 

if(  rpixeKx.y)  ==  WHITE  )  continue; 

if  (  X  ==  0  ) 

{ 

wpixeK  x.y.fid  ); 

continue; 
> 

if(  (  rpixeK  x-l,y  )  !=  WHITE  )  ) 
{ 

wpixeK  x,y,rpixel(x-l,y)  ); 

continue; 
} 

if(  rpixel(x,y-l)  !=  WHITE  ) 
■C 

wpixeK  x.y, rpixeK  x,y-l  )  ); 

continue; 
} 

wpixelC  x.y.fid  ); 
a  =  X  +  1; 

while(l)  .  -    "i 

{ 

if(  rpixeK  a.y  )  ==  WHITE  ) 
{ 

wpixel(  x.y.fid  ); 

step(   ftfid.&gid   ); 

break; 
> 

if(   rpixel(a,y-l)    !=  WHITE   ) 
•C 

wpixelC   x,y, rpixeK   a,y-l   )    ); 

break; 
} 
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a  +=  1; 
} 
}  /*  end  X  -  lor  ♦/ 
>      /»  end  y  -  lor  */ 

lid  +=  gid*HIGH; 

printl("\n\n\tSTART  STEP  TWO...  "); 

printl("\n\nEnter  the  BEST  estimate  ol  the  max  leature  length"); 
geti("\n\nlt's  Better  to  be  too  leirge  than  too  small!  !"  .ftmajcll) ; 
printl("\n\n\tRUHNING") ; 

lor(  y  =  1  ;  y  <  HROW  ;  y++  ) 
{ 

lor(  X  =  1  ;  X  <  NCOL  ;  x++  ) 
{ 

nl  =  rpixeK  x,y  ) ; 

n2  =  rpixelC  x,y-l  ); 

il(  (  nl  I  I  n2  ==  WHITE  )  II  (nl  ==  n2  ))  continue; 

x3  =  X  -  maxll; 
x4  =  X  +  maxll; 
y3  =  y  -  maxll; 
y4  =  y  +  maxll; 

il(  x3  <  LOW)  x3  =  LOW; 
il(  x4  >  NCOL  )  x4  =  NCOL; 
il(  y3  <  LOW)  y3  =  LOW; 
il(  y4  >  NROW  )  y4  =  NROW; 

lor  (  yl  =  y3  ;  yl  <  y4  ;  yl++  ) 
•C 

lor  (  xl  =  x3  ;  xl  <  x4  ;  xl++) 
{ 

nla  =  rpixelC  xl.yl  ); 
il  (  nla  ==  WHITE)  continue; 
il  (  nla  ==  nl  )  wpixeK  xl,yl,n2  ); 
} 
lid  -=  1; 


> 

printl("\n\n\tSTEP  THREE  in  Progress"); 
lor   (   y  =    1    ;    y  <  NROW    ;    y++   ) 

•c 

lor(   X  =    1    ;    X  <  NCOL   ;    x++   ) 
{ 

kk  =  rpixeK   x,y  ) ; 

il(   kk  ==  WHITE  )    continue; 
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kkl  =  kk; 

nnl  =  temp.fid  +  temp_gid  *     HIGH; 

if(   kkl   <  nnl    )    continue; 

if(    (kkl  -  nnl)   <  -350  )  kkl   =  kk  +  508; 

if(    (kkl   -  nnl)   >=   -350     ftft   (kkl   -  nnl)    <  -150   ) 

kkl   =  kk  +  254; 
il(  kkl    !=  nnl   ) 
{ 

x3  =  X  -  majcfl; 

x4  =  X  +  maxfl; 

y4  =   y  +  majcfl; 

if(   x3  <  0)   x3  =  LOW; 

if(   x4  >  HCOL-1    )   x4  =  HCOL  -1; 

if(    y4  >   HROW-1    )    y4   =  HROW  -1; 

for(  yl  =  y   ;   yl  <  y4   ;   yl++   ) 
•C 

for(   xl   =  x3    ;    xl   <  x4    ;    xl++) 
■C 

nlc  =  rpixeK   xl.yl   ); 

if(    (   nlc  ==  WHITE   )    II    (nlc    !=  kk   )    )continue; 
wpixel(   xl.yl ,temp_f id      ); 
} 
> 
step(   &temp_f id,&temp_gid   ); 


printf("\n\n\t\tFEATURE  COUNT  IS:        '/.d'-.fid); 
maplut(    INPUT,    0,    IXS,    lYS,    NCOL,    NROW   ); 
} 
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APPENDIX  G:  PROGRAM  FILE:  sizeit.c 


/♦PURPOSE     This  routine  uses  ain  existing  or  saved  insige  that  has  been 

thresholded  and   then  ID'd  with  the  ID  Algorithm.   The  output 
of  this  program  is  a  tabulair  output  of  the  calculated  Area, 
X-chord,  and  Y-chord  that  is  suitable  for  input  to  a 
Statistical  Analysis  Progrzun.   Output  file  is  written  to 
default  drive  as  SIZE. DAT. 

(  Background  ==  WHITE  ;  Feature  ==  BLACK  ) 

*/ 

#include  "thesis. h" 

char       filename [20] ,  comline[200] ; 
int        srov,  scol; 

mainO 

int  flag; 

printf ("\n\aUSE  IMAGE  processed  by  the  FEATURE  ID  program!!"); 

printf ("\n\tReady  to  Load  ID'd  image  from  disk?... Yes  (y)  /  No  (n)  "); 

flag=getch(); 

if (flag  ==  'Y' I  I  flag  ==  'y') 

{ 

staxtitO  ; 

readitO  ; 
} 

if(  flag  ==   'N'  II  flag  ==  'n'  )  return(O) ; 
initializeO  ; 
size_it() ;  - 

printf ("\n\tSave  image  to  Disk  File?... Yes  (y)  /  No  (n)   "); 

flag=getch() ; 

if  (flag  ==  'Y'll  flag  ==  'y')  saveitO; 

> 

/*  image  feature  sizing  algorithm  */ 

size_it() 
{ 

register  int     nn,  x,  y,  xa,  ya,  xb,  yb; 

int  j,  ip,  ipa,  ipb,  xl ,  yl,  x2,  y2,  maxfl,  num; 
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int 

int 

int 

int 

float 

float 

float 

long 

long 

long 

long 

FILE 


xmax,   ymax,   xmin,   ynin,  maxarea,   mmaxea; 

xflag=0; 

yflag=0; 

aflag=0; 

eirea,  xchord,  ychord,  factor; 

fxmsLX,  fymax,  fxmin,  fymin,  imaxaxea.,   fminEirea; 

cf_l,  cf_2,  mag; 

*aptr; 

*xptr; 

*yptr; 

*stop; 

*f eatdata; 


geti("\nEnter  the  number  of  Features", ftnum) ; 
getf ("\nEnter  Magnification  Factor" .ftmag) ; 
geti("\nEnter  the  Maximum  Feature  Length",fcmaxf 1) ; 

/*  Set  initial  test  min/max  values  */ 


xmeuc  =  0; 

ymax  =  0; 

maxarea  =  0; 

xmin  =  512; 

ymin  =  512; 

minaxea  =  512; 


/*  Dynajnically  allocate  memory  space  for  the  number  of 
features  entered  by  the  operator. 


Note:  This  eliminates  the  need  to  declare  lao^ge  arrays 
for  these  values  eind  allows  maximum  number  of  features 


*/ 


xptr  =  (  long  ♦)  callocC  num  ,sizeof (long)) 
yptr  =  (  long  *)  calloc(  num  .sizeof (long) ) 
aptr  =  (  long  *)  calloc(  num  , sizeof (long)) 


/*  Check  for  successfull  memory  allocation   */ 

if (  !xptr  I  I  !yptr  I  I  !aptr  ) 
•C 

printf ("\nOut  of  Memory!!!  \n"); 

printf ("\nYou  may  have  entered  too  mamy  features  An") ; 

printf ("\nCorrect  problem  and  try  again."); 

exit(O); 
} 

/*  Calculate  Magnification  Constcints        */ 

sccile(&factor) ; 

cf _1   =    (   factor  /  mag  ) ; 

cf_2  =    (   cf_l   *   cf_l    ); 
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/*  Above  Conversion  factor  converts  size  to  microns 

amd  is  related  to  the  input  magnification  for  the 
holograun.   The  defined  "  SCALE_FACTOR  "  constant 
value  must  be  set  for  the  proper  value  based  on  the 
equipment  set  up  during  actued.  image  aquisition. 
The  constant  term  SCALE_FACTOR  is  defined  in  the 
Header  File  THESIS. H 
♦/ 

/*  Begin  Sizing  Routine   */ 

printf ("\n\n\tStarting  to  SIZE!..."); 

nn  =  0;  /*  Feature  number  counter  */ 

for(  y  =  0  ;  y  <  NROW  ;  y++  ) 
{ 

if(  nn  ==  num  )breai;       /*  Quit  when  all  features  sized  */ 

for(  X  =  0  ;  X  <  NCOL  ;  x++  ) 

{ 

ip  =  rpixel(x,y); 

if(  (ip  ==  WHITE)  I  I  (ip  ==  BLACK)  )continue; 

if(  nn  ==  num  )brecik;    /*  Quit  when  all  features  sized  */ 

xl  =  X  -  maxfl;  /*  Set  up  sizing  Box  */ 

yl  =  y; 

x2  =  X  +  maxfl; 

y2  =  y  +  maxfl; 

if(  xl  <  LOW)  xl  =  LOW; 

if(  yl  <  LOW)  yl  =  LOW; 

if(  x2  >  NCOL  )   x2  =  NCOL; 

if(  y2  >  NRGW  )   y2  =  NROW; 

for(  ya  =  yl  ;  ya  <  y2  ;  ya++  ) 
{ 

for(  xa  =  xl  ;  xa  <  x2  ;  xa++  ) 
{ 

ipa  =  rpixel(xa,ya) ; 
if(  ipa  !=  ip  ) continue; 
aflag++; 
xflag++; 
} 

if(  xflag  >  *(xptr  +  nn  )  )  *(xptr  +  nn)  =  xflag; 
xflag  =  0; 
> 

*(aptr  +  nn)  =  aflag; 
aflag  =  0; 
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for(  xb  =  xl  ;  xb  <  x2  ;  xb++) 
•C 

lor(  yb  =  yl  ;  yb  <  y2  ;  yb++  ) 
{ 

ipb  =  rpixel(xb,yb) ; 
il(  ipb  !=  ip  )  continue; 
yflag++; 

wpixeKxb.yb, BLACK) ; 
} 

if(  yllag  >  *(yptr  +  nn  )  )  *(yptr  +  nn)  =  yflag; 
yflag  =  0; 
} 

/*  Calculate  Min/Max  values  */ 

xmax  =  max(  xmax,  *(xptr  +  nn)  ); 

ymax  =  max(  ymax,  *(yptr  +  nn)  ); 

majcaxea  =  meuc(  maxarea,  *(aptr  +  nn)  )  ; 

xmin  =  min(  xmin,  *(xptr  +  nn)  )  ; 

ymin  =  min(  ymin,  *(yptr  +  nn)  )  ; 

minajrea  =  niin(  minairea,  *(aptr  +  nn)  )  ; 

nn++;  /*  increment  counter  */ 

printl("\n\n\tFeature  '/.d  complete.  .."  ,nn)  ; 
} 
}      /*....  end  initial  x/y  for  loop  ....*/ 

/*  output  data  section    */ 

printf("\n\a  Sending  output  data  to  screen  emd  FILE  SIZE. DAT  !"); 
/*  waitO;   */ 

if(  (  featdata  =  f open("size.dat" ,"w")  )  ==  NULL  ) 
{ 

printf  ("CANNOT  Open  output  file  '/.sXn",  *f  eatdata)  ; 

return(l) ; 
} 

printf ("     ID  NO       AREA    X-Width   Y-Width  \n"); 

for(  j  =  0  ;  j  <  num  ;  j++  ) 
{ 

fprintfC  featdata, '"/.lOd  '/,10.3f  7,10. 3f  '/.lO.SfXn",  j  +  1,  *(aptr  +  j)  *  cf_2, 

*(xptr  +  j)  *  cf_l,  *(yptr  +  j)  *  cf_l  ); 

printf ("\n'/.10d  '/.lO.Sf  '/.lO.Sf  '/.lO.Sf  ".  j  +  1,  *(aptr  +  j)  *  cf_2, 
*(xptr  +  j)  *  cf_l,  *(yptr  +  j)  *  cf_l  ); 
} 

/*  This  section  csin  be  used  to  put  total  information  on  bottom  of  data  file 
if  desired.   File  written  now  to  be  used  with  STATGRAPHICS  */ 
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/*  lprintf(featdata,"\n\n\nMax  X-chord=  '/.f  Max  Y-chord=  V.f  Max  Area=  '/.f   "\ 

xmax*ct_l  ,yinaLX*cf_l  ,maxajrea*cf  _2)  ; 

fprintf(featdata,"\nMin  X-chord=  V.i         Min  Y-chord=  '/,i         Min  Area=  V.f  ", 
xinin*cl_l,3rmin*cf_l,minaxea*cf  _2) ;   */ 

printl("\n\n\iiMax  X-chord=  '/.i         Max  Y-chord=  */,!         Max  Area=  '/,f   ", 
xmax*ci_l  ,yniaLX*cl_l  ,m2Lxajrea*cl_2) ; 

printfC'XnMin  X-chord=  '/.i         Min  Y-chord=  V.f         Min  Area=  '/.f   ", 
xinin*cf  _1 ,3rmin*cf  _1  ,minzurea*cf  _2) ; 

/*  Close  open  file  and   free  allocated  memory   */ 

f close(leatdata) ; 
Iree(xptr) ; 
Iree(yptr) ; 
free(aptr); 


62 


APPENDIX  H:  PROGRAM  FILE:  speckle.c 


/*PURPOSE   :  Provides  for  calculating  the  speckle  index  of  an   image. 
The  resulting  value  is  used  as  a  measure  to  evaluate  the 
sucess  of  filtering  an   image  to  reduce  the  speckle  noise. 
This  routine  is  called  in  all  filter  eilgorithms . 

— >  OPERATOR  must  enter  a  lumber  of  ROWS  and  COLUMNS 
THE  PROGRAM  uses  to  calclate  the  speckle  index. 

ROWS  >  Maximum  =  480      Minimum  =  1 

COLS  >  Maximum  =  512      Minimum  =  1 


*/ 

#include  "thesis. h" 

speckleO 

{ 

extern  int 

srow,  scol; 

int 

m,  n,  X,  y, 

nn; 

long 

smax,  smin; 

long 

sdataCNUM] ; 

unsigned  long  deviation,  ssum; 

float  smeam,  slocal,  stotal,  tot; 
float  spklindex; 

if(  sro¥  >=  NROW  )  srow  =  NROW; 
if(  scol  >=  NCOL  )  scol  =  NCOL; 

tot  =  (long)srow  *  scol; 

/*  Commence  calculation   */ 

printf ("\n\n\tCalculating  Speckle  INDEX..."); 

for  (  n  =  1  ;  n  <  srow-2  ;  n++) 
{ 

for  (  m  =  1  ;  m  <  scol-2  ;  m++) 
{ 

smax  =  0 ; 
smin  =  260; 
ssum  =  0; 
nn  =  0; 
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lor(  y  =  (  n-1  )  ;  y  <  (n+2)  ;  y++  ) 
{ 

lor(  X  =  (  m-1  )  ;  X  <  (  m+2  )  ;  x++) 
{ 

sdataCnn]  =  rpixel(  x,y  ); 
ssum  +=  (long)sdata[nn] ; 
if  (  smeix  <  sdata[im]  )  smajc  =  sdataCnii]  ; 
il(  smin  >  sdata[im]  )  smin  =  sdataCnn]  ; 
nn++; 
} 
} 

deviation  =  smax  -  smin; 
if(  ssum  ==  0  )  smeem  =  1; 
smean  =  ssum  /  HUM   ; 
slocal  =  (lloat)deviation  /  smeeoi; 
stotal  +=  slocal; 
} 
} 

spklindex  =  stotail  /  tot; 

printf  ("\n\n\tThe  Calculated  */,d  by  V.d  speckle  index  is  */,f  "  .srow.scol, spklindex) ; 
} 

/* END  Speckle  Function  */ 
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APPENDIX  I:  PROGRAM  FILE:  vir_arry.c 


/♦PURPOSE   :  These  functions  provide  lor  the  setup  ol  a  disk  drive  virtual 
cirray  that  can  be  indexed  as  if  it  vere  in  the  calling 
programs  data  storeige  axea.   The  routines  provide  for  disk 
file  access  when  required  to  retrieve  data  elements  from  the 
array.  This  allows  the  size  of  a  declared  array  to  be  limited 
only  by  the  amount  of  DISK  SPACE  available.  See  "thesis. h" 
header  file  for  setup  and   HACRO  routines  that  support  these 
functions.   Also  progrsim  Istat.c  fully  implements  the  eurray 
routines . 

♦/ 
#include  "thesis. h" 

/*  Virtual  Array  Access  Routines  */ 

int  init_v_2irray(f  ilename,rec_size,f  ilchax) 

chcLT  *filename,    filchax; 

int  rec_size; 


long    size; 

FILE    *f; 

f  =  fopen(f  ilensone/'wb")  ; 

if  (f  !=  NULL) 

{ 

size  =  0; 

fwrite(&size,4,l,f ) ;  /*  write  array  size  of  zero  */ 
fwrite(&rec_size,2,l,f ) ;  /*  and  array  element  size  */ 
fwrite(&filchar,l,l,f);      /*  and  fill  char  */ 

fclose(f);  /*  to  file  header  ♦/ 

return(l) ; 


} 

else 


return(NULL); 


VACB  *open_v_aLrray  (filename,  buff  er_size) 

chcir  *filenaine; 

int  buffer_size; 
{ 

VACB    *v_array; 

char    *buf_ptr; 

int     i; 

chax    filchar; 
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/*  allocate  virtual  array  control  block  */ 

v.array  =    (VACB  *)   malloc(sizeof (VACB)) ; 
if   (v_array  ==  MULL)   return(HULL) ; 

/*  open  virtual  array  file  */ 

v_array->f ile  =  fopen(f ileneUDe,"r+b") ; 

if    (v_array->file  ==  HULL) 

{ 

free(v_aurray) ; 

return (MULL); 

}; 

/*  get   eorray  size  aind  element   size  for  control  block  */ 

fread(ftv_eirray->size,4, 1  ,v_ea:ray->f  ile) ; 
fread(&v_axray->elsize,2, 1 ,v_2a-ray->f ile) ; 
f  read(ftf  ilchax  ,1,1,  v_cLrray->f  ile)  ; 
v_aarray->buf_elsize  =  v_axray->elsize  +  4; 

/*  allocate  buffer  ♦/  , 


v_array->buff er  =   (char  *)   malloc(v_2u:ray->buf_elsize   *    (buffer_size   +   1)); 

if  (v_array->buffer  ==  NULL) 

{ 

f  close(v_au:ray->f  ile) ; 

f ree(v_array) ; 

return (NULL); 

}; 

v_array->buf_size   =  buffer_size; 

/*  set  up  blank_rec  using  the  fill   character   in  eirray  header  */ 
/*  for   initializing  new  ao-ray   elements   */ 

buf_ptr  =  v_aLrray->buff er  +  v_array->buf _elsize  *  v_aLrray->buf_size; 
v_cirray->blcink_rec  =  buf_ptr  +  4; 
for   (i   =  0;    i  <  v_axray->buf _elsize;    i++) 
*buf_ptr++  =  filchar; 

/*  set  record  index  negative  for  all  buffer  elements  */ 

buf_ptr  =  v_array->buff er ; 

for   (i   =  0;    i   <  v_eLrray->buf _size;    i++) 

{ 

*((long  *)buf_ptr)    =   -IL; 

buf_ptr  +=  v_array->buf _elsize; 

>; 

return  (v_cirray) ; 
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void 
VACB 


close_v_array(v_axray) 
*v_aj:ray; 

int         i; 

cheoc       ♦bul.ptr; 

long       rec_index,   lile.ollset; 

bul_ptr  =  v_aurray->buff  er; 


/*  Hush  buffer  ♦/ 

for   (i  =  0;    i  <  v_axray->buf_8ize;    i++) 


/*  check  each  element  index  */ 

/*  if  element  present;  write  it  to  disk  */ 

rec_index  =  *((long  *)buf_ptr); 

if  (rec_index  >=  0) 

{ 

file_offset  =  header  +  rec_index  *  v_aLrray->elsize; 

f  seek (v_eirray->file,f  ile_off  set  ,0) ; 

f write (buf_ptr  +  4,    v_axray->elsize,l ,    v_axray->f ile) ; 

}; 

buf_ptr  +=  v_array->buf _elsize; 


}; 

free(v_array->buff er) ; 
f close(v_axray->f ile) ; 
free(v_array) ; 


/*  de-allocate  buffer  */ 
/*  close  array  file  */ 
/*  de-aaiocate  VACB       */ 


void 

VACB 

long  index ; 

{ 


*access_v_rec( v_axray , index) 
*v_array; 


char    *buf_ptr; 

int     buf_index; 

long  rec_index,  temp_index; 


/*  calculate  buffer  address  of  referenced  element  */ 

buf_ index  =  index  '/,  v_array->buf_size; 

buf_ptr  =  v_aj:ray->buff er  +  buf .index  *  v_array->buf_elsize; 

rec_index  =  *(long  *)buf_ptr; 
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/*  it   element  present,  return  its  buffer  address  */ 
a   (rec_index  ==  index)  return(buf _ptr  +  4); 

/*  if  element  doesn't  exist,  extend  the  file  */ 

if   (index  >=  v_eurray->size) 
{ 

f seek (v_aTray->file, 0,2) ; 

for  (temp_ index  =  v_array->size;  temp_index++  <=  index;  ) 

fwrite(v_2Lrray->bl2aik_rec,   v_array->elsize,    1,   v_zun:ay->f  ile) ; 

v_aarray->size  =   index  +   1; 

f seek (v_2Lrray->file, 0,0)  ; 

fHrite(*v_axray->size,  4,  1,  v_array->f ile) ; 

}; 

/*  if  buffer  slot  is  occupied  by  another  element,  ♦/ 
/*  save  it  to  disk   */ 

if  (rec_index  >=  0) 
{ 

f seek(v_2Lrray->f ile,    rec_index   *  v_array->elsize  +  header,    0); 

f writeCbuf _ptr  +  4,    v_axray->elsize,    1,    v_aj:ray->f ile) ; 

}; 

/*  read  referenced  element  into  buffer  slot  */ 

f seek(v_axray->f ile,  index  ♦  v_array->elsize  +  header,  0); 
fread(buf _ptr  +  4,  v_2a:ray->elsize,  1,  v_axray->f ile) ; 
*((long  *)buf_ptr)  =  index; 

/*  return  address  of  element  in  buffer  */ 

return(buf _ptr  +  4) ; 
> 

/♦ */ 
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APPENDIX  J:  PROGRAM  FILE:  genfunc.c 


/♦PURPOSE   :  These  general  fxxnctions  provide  for  routine  evolutions 
that  occur  a  number  of  times  in  the  Image  processing 
routines  developed  for  analysis  of  speckle  reduction 
algorithms.  Some  require  support  of  other  functions 
and  sire  not  totally  independant . 

*/ 
#include  "thesis. h" 
/*  initial  ITEXPC  Board  setup  */ 

startitO 
{ 

sethdw(REGBASE,MEMBASE, MEMORY) ; 

setdim(XSIZE,YSIZE, DEPTH) ; 

aclear (IXS . lYS , NCOL , NROW , COLOR) ; 

fgonO; 

initialize  0 ; 

clsO; 

select_m6m(MEM_A) ; 

display _mem(MEM_A) ; 


/*  Get  integer  keyboard  input   */ 

geti(name,iptr) 
char   *name; 
int   *iptr; 
{ 

printfC"  '/,s  ",name); 

scanfC  '/.d",iptr); 
} 


/*  Get  floating  point  keyboard  input   ♦/ 

getf (rname.riptr) 
chcir   ♦mame; 
int   *riptr; 
{ 

printf("'/,s  ".rname); 

sceinfC"  '/,f  "  .riptr) ; 
} 


/*  Pause  routine  used  to  prompt  for  operator  intervention*/ 
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HaitO 
{ 

printf("\n\n  Press   Return  to  continue"); 

fflush(stdin) ; 

getchJurO; 
} 


/*  returns  minimum  integer  value  from  input  airray  passed  */ 

getmin(data,num) 

int   *data; 
int   num ; 
{ 

int  t,  min; 

for    (min=   dataCO],    t=l    ;    t<nuin    ;    t++) 

if(   data[t]    <  min  )   min  =  data[t] ; 

return  min; 
} 

/*  returns  majcimum  integer  value  from  input  axray  passed  */ 

getmax (data , num) 

int  *data; 
int  num; 
{ 

int  t ,  max ; 

for  (max=  data[0] ,  t=l  ;  t<num  ;  t++) 

if(  dataCt]  >  max  )  max  =  data[t]  ; 

return  max; 
}  - 

/*  DOS  System  Call  to  clear  Display  Screen  */ 

clsO 
{ 

_clearscreen(_GCLEARSCREEN) ; 
} 

/*  routine  to  read  an  image  from  disk  */ 

readitO 

/♦program  uses  the  following  externals:  filename , comline*/ 

{ 

extern  char   filename [20] ,  comline[200] ; 


int  errval ; 

int         ch; 

while(l) 
{ 

printfC'XnXnENTER  IMAGE  FILEHAME  as  (DEV : name. IMG) ->") ; 

lllushCstdin) ; 

scamf  ("'/,s",lilenciine) ; 

errval   =  readiin(IXS,IYS,NCOL,HROW,lilenajne,comline) ; 

if (errval  ==  0)  { 

printf  ("\nFILE->  '/.s  \n\nCOMMEHT :  '/.s"  .filename,  condine) ; 
breaik; 
} 

if (errval  !=  0)  { 
SB itch (errval) 
{ 

case  FILE_ERROR: 

printf ("Error  opening  file\n"); 
breeik; 
case  FORMAT.ERROR : 

printf ("Unknown  file  format\n"); 
breai; 
case  READ_ERROR: 

printf ("Error  Reading  file\n"); 
break; 
default : 

printf  ("?  Unknown  Error  ?'/,d\n"  .errval)  ; 
break; 
} 

printf ("Try  Again??  . . .Yes(y)  /  No(n)\n"); 
ch  =  getchO  ; 

if(ch  ==  'Y'  II  ch  ==  'y')  continue; 
if(ch  ==  'N'  II  ch  ==  'n')  exit(O); 


/*  routine  to  save  an  image  to  disk  ♦/ 

saveitO 

/♦program  uses  the  following  externals:  *f ilename.*comline*/ 

{ 

extern  char       f  ilencime[20]  ,    comline[200]  ; 

int  format,    errval; 

int  ch; 

initializeO ; 

while(l) 
{ 
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printlC'\n\nEHTER  FILENAME. . . (DEV:name. IMG)->") ; 

ff lush(stdin) ; 

sceinf  ("'/,s"  .filename) ; 

geti("\n  Enter  File  Compression  Value  ?. . . (0/l)",*forTnat) ; 

printf("\n\n  Enter  Image  Comment  ..to  continue..  (Max  200  CHAR)..."); 

Iflush(stdin) ; 

gets(comline) ; 

if (format  ==  l)printf ("\n\nStoring  Imaige  Using  Compression  —  Please  wait  !!"); 
if (format  ==  0)printf ("\nStoring  Image  —  Please  wait!!"); 

errveil   =  saveiffl(IXS, IYS,NCOL,NROW, format  .fileneune, comline)  ; 

if(errval  ==   0)      { 

printf ("\nlmage  save  completed  Sat"); 

breaik; 
} 

if(errval    !=  0)      { 

printf ("\n\a  Error  saving  file!!"); 
if(errval  ==  ALLOCATION_ERROR)  { 

printf ("\nlnsufficient  Disk  Space"); 
} 
if(errval  ==  WRITE.ERROR)  { 

printf ("\nError  writing  file  or  values"); 
} 

printf ("\n\tTry  Again??  ...Yes(y)  /  No(n)"); 
ch  =  getchO  ; 

if(ch  ==  'Y'  II  ch  ==  'y')  continue; 
if(ch  ==  'N'  II  ch  ==  'n')  break; 
} 
} 
} 

/*  counter  step  routine  */ 

step(u, v) 
int  *u ,  *v ; 
{ 

*u  +=  1; 

if(*u  >=  HIGH)   { 
*u  =  1; 
*v  +=  1; 
> 
} 

/*  threshold  routine  for  feat_id.c*/ 

threshO 
{ 

int  aiis  ,c,option; 
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unsigned  int  val; 

Bhile(l) 
{ 

initlutsO ; 

initializeO ; 

geti("\n\n\tEirrER  HEW  THRESHOLD  LEVEL  (0-2SS)  :",ftval); 

setlutClHPUT.O); 

thresholdC INPUT, 0, HIGHEST. val) ; 

setlut(RED,0); 

thresholdCRED.O, HIGHEST, val) ; 

8etlut(GREEH,0); 

threshold(GREEN , 0 , HIGHEST , val) ; 

setlut(BLUE,0): 

thresholdCBLUE , 0 , HIGHEST , val) ; 

setlut ( INPUT, 0); 

printf  ("\n\nSATISFIED  WITH  THRESHOLD  AT  '/.d?.  .  .Yes(y)  /  No(n)  .  .  . "  ,val) ; 

ans  =  getchO  ; 

if(ans  ==  'Y'  II  ans  ==  'y')  break; 
} 

maplut (INPUT, 0,IXS,IYS,NCOL,NROW); 
initializeO ; 
} 

/*  threshold  routine  for  threshit.c  */ 

threshitO 
{ 

int   2ais,c,option; 

int   *val ; 

while(l) 
{ 

initlutsO  ; 

initializeO; 

geti("\n\nENTER  NEW  THRESHOLD  LEVEL  (0-255)  :",ftval); 

printf("'/,d\n",val); 

setlut ( INPUT. 0); 

threshold( INPUT. 0. HIGHEST, val)  ; 

setlut (RED, 0);  "  ' 

threshold(RED,0, HIGHEST, val) ; 

setlut (GREEN, 0); 

threshold(GREEN , 0 , HIGHEST . val) ; 

setlut (BLUE, 0);  ,  ' 

threshold(BLUE.O. HIGHEST. val); 

setlut ( INPUT. 0); 

printf  ("\n\nSATISFIED  WITH  THRESHOLD  AT  '/.d?  .  .  .  Yes(y)  /  No(n)  .  .  .  "  ,  val) ; 
cms  =  getchO  ; 
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il(ans  ==  'Y'  II  ans  ==  'y')  break; 
} 

printf("\n\t\t  0:  SAVE  THE  THRESHOLD  IMAGE?..."); 

printl("\n\t\t  1:  QUIT  (Leave  thresholded  image  lor  further  processing)..."); 
printl("\n\t\t  2:  RESTORE  system  to  original  input  image..."); 
geti("\n\t  Select  option  by  lUMBER:  ".ftoption); 
8Bitch(option) 
{ 

case  0: 

maplut(INPUT,0,IXS,IYS,ICOL,IROW); 
saveitO ; 
bresik; 
case  1 : 

maplut ( IIPUT , 0 . IXS , lYS , SCOL , MROW) ; 
breeiX; 
case  2: 

initializeO ;       /*initialize  without  mapping  operation*/ 
break; 
default : 
breeik; 
> 
} 

/♦   function  to  compute  standard  deviation  of  image   */ 

dev(stddev) 

float       *stddev; 
{ 

float       var,  sum,  sqsum,  sqvalue,  dim,  value; 

long       X,  y; 

sum  =  sqsum  =  0; 

dim  =  (float)NROW  *  NCOL; 

for(  y  =  0  ;  y  <=  NROW-1  ;  y++) 

{ 

for(  X  =  0  ;  X  <=  NCOL-1  ;  x++) 
{ 

value  =  rpixel(x,y); 
sqvalue  =  value  *  value; 
sum  +=  value/dim; 
sqsum  +=  sqvalue/dim; 
> 
} 

var  =  sqsum  -  (sum*suin); 
*stddev  =  sqrt(var)  ; 


/*  function  to  compute  scale  factor  for  sizing  routine   */ 
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sc2Lle(lactor) 

float        *1 actor; 


int  inch,  pixel,  flag; 

float  dinch,  dpixel;      /*  default  length  values  */ 

/*  dinch  amd  dpixel  aire  to  be  modified  accordingly  ♦/ 

dinch  =   1 . 00 ; 

dpixel  =  50.00; 

printf("\n\tUtili2ing  standard  SCALE.FACTOR?. . .Yes  (y)  /  Ho  (n)"); 

flag  =  getchO; 

if(flag  ==  'Y'  I  I  flag  ==  'y') 

{ 

♦factor  =  ((float)SCALE_FACTOR); 
} 

if (flag  ==  'N'  II  flag  ==  'n') 
{ 

printf ("\n\tPress  '0'  if  using  default  values  for  length  in  inches  or  pixels"); 

geti("\n\tEnter  the  length  of  object  in  inches(def ault) : ".ftinch) ; 

geti("\n\tEnter  the  length  of  object  in  pixels(def ault) : ".ftpixel)  ; 

if (inch  ==  0) 
{ 

inch  =  dinch; 
} 

if(pixel  ==  0) 
{ 

pixel  =  dpixel; 
} 
♦factor  =  ((float)inch)  /  ((float)pixel) ; 

> 

printf  ( "\n\tSCALE_FACTOR  is  y.f\n" ,  ♦factor)  ; 
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APPENDIX  K:  C  COMPILER  OPERATIONS 

The  following  batch  files  Hst  the  C  compiler  directives  used  to  compile  and  link 
all  C  programs.  The  C  compiler  files  ACCOMP.BAT  and  DACCOMP.BAT  are 
listed  below.  ACCOMP.BAT  contains  the  switches  necessary  to  utilize  CODEVIEW 
debugger,  and  ACC0MP2.BAT  is  for  final  program  optimization.  The  F  drive  Hsted 
in  the  PATH  statement  of  these  batch  files  is  due  to  the  fact  that  the  programs  were 
stored  on  a  Bernoulli  disk. 


ACCOMP.BAT: 

path  c:\; c:\dos; c:\wsjet ;f :\kedit .*;c:\mousel ;c:\xtree; c:\lsrctrl; 

1 :\;f :\msc\bin;f :\msc\tmp;f :\msc\src;f :\msc\wrk;f :\msc\lib; . . 

l:\msc\iiiclude 
set  include=\msc\include 
set  lib=\msc\lib;\pcplus\itex\lib 
set  tmp=\msc\tmp 

cl  /c  /Fs  /Zi  /GO  /AL  /Fm  '/.I .  c 

ACC0MP2.BAT 

path  c:\; c:\dos; c:\wsjet ;f :\kedit .*;c:\mousel ;c:\xtree;c:\lsrctrl; 

l:\;f :\msc\bin;f :\msc\tmp;f :\insc\src;f :\msc\wrk;f :\msc\lib; . . 

f:\msc\include 
set  include=\msc\include 
set  lib=\msc\lib;\pcplus\itex\lib 
set  tmp=\msc\tmp 

cl  /c  /Fs  /Ze  /GO  /AL  '/.l.c 


The  following  batch  files,  LTHESIS.BAT  and  DTHESIS.BAT,  Hnk  the  programs 
to  the  libraries  and  supporting  function  files.  The  speckle. c  should  not  be  included 
within  these  files  when  linking  the  programs  sizeit.c  and  feat.id.c  due  to  an  undefined 
external  problem  associated  with  the  external  variables  SROW  and  SCOL.  LTHE- 
SIS.BAT is  used  for  debugging  applications,  and  DTHESIS.BAT  is  for  final  program 
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optimization.  It  is  also  unnecessary  to  link  feat-id. c  and  sizeit.c  to  the  virtual  array 
functions  contained  within  vir^arry.c. 


LTHESIS.BAT: 

IF  HOT  ERRORLEVEL  1  LINK  /lOD  /CO  /E  '/.l  genlunc   speckle  vir.arry 
, , , ITEXPCML . LIB  LLIBCE . LIB ; 


DTHESIS.BAT: 

IF  HOT  ERRORLEVEL  1  LINK  /HOD  /E  '/.l  genfunc   speckle  vir.arry, 
, , , ITEXPCML . LIB  LLIBCE . LIB ; 


To  create  the  menu  format  for  program  display,  a  C  program  named  hologram. c 

was  generated.    To  link  all  image  processing  programs  to  this  module,  all  program 

"main()"  statements  were  modified  so  that  hologram. c  was  made  the  main  program. 

All  image  processing  programs  were  then  hnked  to  hologram. c  after  compilation.  The 

ACC0MP2.BAT  file  was  used  for  compilation,  and  the  PTHESIS.BAT  file  was  used 

for  subsequent  module  linking.    Due  to  the  length  of  the  LINK  statement,  the  user 

must  enter  the  libraries  ITEXPCML  and  LLIBCE  manually  into  the  keyboard  each 

time  the  LINK  process  is  accomplished. 

PTHESIS.BAT: 

IF  NOT  ERRORLEVEL  1  LINK  /NOD  /E  '/.l  table  speckle  vir.arry... 

genfiinc  threshit  sizeit  feat.id  geof  il.  . 

Istat  2sigina 
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